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The Effect of Jet Noise on Aircraft Structures 


B. L. CLARKSON, B.Sc., Ph.D. 


(Department of Aeronautical Engineering, The University of Southampton) 


SumMARY: The present state of knowledge on the problem of fatigue failure due 

to vibrations excited by jet noise is reviewed. It is concluded that it should 

currently be possible to make reasonable estimates of the stress levels set up 

in a structure by jet noise but, in general, the resultant fatigue life of the 
components cannot be estimated with any confidence. 


1. Introduction 

The efflux from a jet engine produces fluctuating noise pressures which cause 
vibration of those parts of the aircraft structure situated in the region of the jet. 
The pressures are varying randomly in the statistical sense and cover a wide range 
of frequencies. Therefore, as a result of this wide frequency spectrum, resonant 
vibrations are set up in the structure. Such vibrations are generally of low 
amplitude but, because of their high frequency, a large number of stress reversals 
can be produced in a relatively short running time. The possibility of fatigue 
failure of these parts must now be considered in the design of the structure. The 
overall solution of the problem requires a study of many interrelated topics and the 
purpose of this paper is to review and co-ordinate the work being done in this field. 


The complete analysis of the problem of fatigue failure due to vibrations 
excited by jet noise requires a detailed study of each of the following items :— 


(i) General theory of structural response to random pressure loads. 

(ii) Characteristics of the noise pressure field. 

(iii) Mechanical impedance of the aircraft structure. 

(iv) Practical approximations to the general case by considering panel 

vibrations. 

(v) Fatigue life estimation. 

The applied loads on the structure take the form of randomly varying pressures 
associated with noise which has a wide frequency spectrum. A statistical approach 
to the problem by using the method of generalised harmonic analysis is therefore 
the most convenient. By considering each normal mode of vibration of the 


structure separately, it is possible to use Lagrange’s equations of motion to obtain 
the mean square value of the response to random loads. The random loading 
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must be defined by its power spectrum and also by its spatial correlation over the 
structure. The effectiveness of this loading in exciting a given mode of vibration 
is a function of its spatial correlation and of the mode deflection. The response of 
a structure to a given loading depends on the vibration characteristics of the 
structure. In aircraft structures the damping is low and the power spectrum of the 
noise pressure is relatively flat. Thus a practical solution can be obtained on the 
assumption that the response is essentially resonant. 


It has been shown that the relevant properties of the noise pressure field are 
the pressure power spectrum and the spatial correlation of pressures over the whole 
structure. The mechanism of noise production by a jet efflux is so complicated 
that these properties must be obtained experimentally. The power spectrum of jet 
noise is relatively flat, indicating that the random pressure fluctuation can be con- 
sidered to be made up of components in a wide range of frequencies. Close to the 
nozzle the power spectrum peaks towards the higher frequencies, while farther 
downstream the peak moves towards the lower frequencies. There is as yet no 
means of calculating this power spectrum to within the accuracy required for 
structural response calculations. Because of random variations in the noise 
pressure field close to the jet, the phase relationship along the structure must be 
obtained by correlation measurements. In this region the area of positive correla- 
tion of the overall pressure is of the order of one nozzle diameter. For points in 
the far field the jet can be considered to form a point source and hence the 
correlation can be calculated directly. 


The vibration characteristics of the structure are defined by the normal modes 
and natural frequencies of vibration, together with the associated damping. In the 
present state of knowledge the normal modes of vibration of the general case of a 
stiffened cylinder of the type used in aircraft construction cannot be determined 
theoretically. However, in the special case of stiffeners, whose spacing is close 
compared with the flexural wavelengths, a modification of the simple case of a 
uniform cylinder can be made to render the problem tractable. This assumption 
of closely spaced longitudinal and circumferential stiffeners may adequately cover 
about half the number of modes significant in jet noise excitation. Thus further 
investigation is urgently required into the general case. The damping of the normal 
modes is due to energy dissipation in the material itself, joint slip, and sound 
radiation. As the mean square stress level is inversely proportional to the damp- 
ing ratio, an increase in the damping will automatically produce a reduction in 
stress levels. Structural damping can be increased by the addition of damping 
materials to the inner surface of the structure or as interlayers at the joints. 
Acoustic damping cannot be varied appreciably on a given aircraft, but in the 
design stage of a new aircraft it may be possible to design for high acoustic 
damping. 


As the correlation area is of the same order as the area of typical panels, it is 
probable that the form of vibration of the complete structure approximates to that 
of predominantly panel or stringer vibrations. However, this is a special case of 
the overall vibration of a stiffened cylinder subjected to localised loading and should 
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EFFECT OF NOISE ON STRUCTURES 


be treated generally until experience on typical structures justifies the approxima- 
tion. Considerable experimental and theoretical work has been done on single 
panels subjected to high intensity acoustic loading, but in general it is not possible 
to extrapolate these results to the case of the complete structure. Such tests provide 
an indication of the relative fatigue lives of different panels but cannot reliably 
predict the magnitude of the fatigue life on the actual aircraft. 


It can be concluded, therefore, that in the present state of knowledge it should 
be possible to make reasonable estimates of the stress levels set up in a structure 
by jet noise, but the fatigue life of components under these conditions cannot be 
estimated with any confidence. The effect of random loading at low stress levels 
needs considerable investigation. However, in special cases it may be possible to 
make estimates based on full-scale test experience on a similar structure and then 
use the theory to allow for the effect of structural changes. 


NOTATION 


area of structure, amplitude of axial displacement 


element of area 


A 
a 
B amplitude of circumferential displacement 
b panel wavelength 
C amplitude of radial displacement 
c generalised viscous damping, velocity of sound 
f. equivalent frequency (c./s.) 
h_ constant to convert displacement to stress 
j joint acceptance 
k_ generalised stiffness 
L generalised force 
I length of cylinder 
m_ generalised mass, number of longitudinal half-waves 
N’ number of root-mean-square stress reversals to produce failure 
total number of modes, number of circumferential waves 
overall noise pressure 


noise pressure in a narrow frequency band 


n 
Pp 

R_ correlation coefficient 

S panel stress due to static pressure 
T time to failure 

t time 


u_ axial displacement 
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circumferential displacement 

radial displacement 

distance along cylinder 

mechanical impedance 

mode deflection, ratio of stress to root-mean-square stress 
damping ratio 

reduction factor 

angle of incident sound waves 

sound wavelength 

generalised co-ordinate 


ore Ne =F 


Stress 


~ 


time delay 


radial position 


€ 


frequency (rad./sec.) 


Suffixes 
0 denotes reference point 


Ss denotes general mode 


2. General Theory of Structural Vibrations Excited by Random 
Loads 


2.1. INTRODUCTION 


When a structure is subjected to random pressure loads, associated with a jet 
efflux, the resultant stresses are also random in nature. In this case instantaneous 
values of stress are meaningless, and therefore the problem must be treated 
Statistically. The pressure loading has a wide frequency spectrum, causing excita- 
tion of many natural modes of vibration of the structure. In these circumstances 
the most convenient method of analysis is to use the Lagrange equations to consider 
the excitation of each normal mode of the structure independently. This assumes 
that the coupling between modes is negligible, and is justified by the fact that the 
damping of the normal modes is small. The total response of the structure is then 
obtained by combining the effects in each mode. The general study of random 
loading of a structure indicates what information about the exciting forces is 
required and shows in general terms the form of the response. It is shown later 
that, owing to the peculiarities of the system and of the excitation, it is possible to 
make simplifications which render the problem tractable. The method used is that 
of Generalised Harmonic Analysis developed by Liepmann® and others and 
applied by Powell’: *: *’ to this particular problem. 
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2.2. RESPONSE OF A STRUCTURE TO RANDOM LOADING 


Consider the response of the structure, in the s“* mode of vibration due to 
random noise pressures. The Lagrange equation for this mode has the form 


Meet (s=1 tu n) 


where m,, C,, ks, L,(t) are the generalised mass, generalised viscous damping, 
generalised stiffness and generalised force respectively. £€ is the generalised 
co-ordinate of the mode, i.e. the deflection of the reference point. 


The solution of this equation yields the power spectral density of the 
displacement for the s“ mode, 


Wz (w) 


We (o)= Y |?’ 


and the mean square value of displacement is 
= Wz (@) 
We (w) do= ¥ do, 


0 


where Y () is the mechanical impedance of the system, i.e. 
Cc? 
[¥ P=m | 


,=natural frequency of the s mode 


and wy, (w#)=power spectral density of the generalised force. 


Now the generalised force 
| P(.a)y(@ da, 
A 


where P(t, a) is the randomly varying noise pressure on the element of area of the 
structure, and y(a) is the mode deflection relative to the reference point, which is 
taken as unity. The square of the generalised force is 


|| (a) dada. 


AA 


In this expression the product P(t, a) P(t, a’) is the spatial correlation of the total 
noise pressures at positions a and a’. 
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The power spectral density of the generalised force at frequency » is 


w.(o)= | | papa) y@y (a) dada, 


AA 


where p (t, a) p (t, a’) is now the spatial correlation of the noise pressures in a narrow 
frequency band of central frequency ». This correlation function can be normalised 
by dividing by the mean square pressure at some reference element ay. 


Re-writing : — 


p(t, a) p(t, 


AA 


where now 


p(t, a) p(t,@)=the correlation coefficient in a narrow band having a 


P(t, a) central frequency 
=R., say 


and p(t, a,)*=the power spectral density of the noise pressures in a 
narrow band at frequency » and station a, 


=Wo (o), say 


Therefore w, (o)= | Wyo y (a) y da da’. 


AA 
Thus the power spectral density of displacement 


wif) 
| Woo () Ruy (a) y (a) dada’. 


AA 


We 


As the mode shape is known, the power spectral density of the stress at any 
point can be determined directly from w;(w). It is simply 


We (w)= 2 j Wyo Ruy (a) y (a) da dd, 
A 


A 


where h, is a constant converting the displacement of the reference point in the s* 
mode to the stress at any required point in the mode. 
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It is seen from these expressions that the information required to arrive at a 
solution is 
(i) Wp (w), the power spectral density of the noise pressure at some reference 
point. 
(ii) R., the narrow band spatial correlation of the noise pressures over the 
whole structure. 


(iii) Y (), the mechanical impedance of the structure which is dependent on 


(a) The normal structural modes and natural frequencies of vibration. 
(b) The damping of the normal modes. 


(iv) y(a), the mode shape. 


2.3. FURTHER CONSIDERATION OF THE GENERALISED FORCE 
2.3.1. General 

Several specialised cases of the generalised force can now be considered. It is 
convenient to re-write the expression for the generalised force as 


W, (w) A? | R.y (@ y da dd | 
AA 


The term in square brackets is now non-dimensional, having values equal to, or 
less than, unity. It is a measure of the inter-relationship between pressures and the 
mode deflections and its square root has been termed “ joint acceptance ” by Powell 
and denoted by j. We have then simply that the power spectral density of the 
generalised force is given by 


WL (o)= Woo () 


where j is a factor to be determined from the space correlation and the mode 
deflection. 


2.3.2. Near Field 

For structure in the near field of the jet it is not possible to calculate R.. 
Thus j? must be obtained from experimentally determined values of R.. In the 
special case of a correlation coefficient falling off rapidly to zero in a distance which 
is short compared with the half wavelength of the mode, the mode deflections 
y (a), y (a) will be approximately equal over the region having appreciable values 
of R.. 


2 
Then 7? (0) R.y? (a) da= ~ R.da, 
A A 
which can readily be evaluated. When this condition is not true the calculation of 
? (@) is best carried out on a digital computer. 
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te 


FiGcureE 1. Pressure loads on a panel due to an incident sound wave. 


2.3.3. Far Field 


For structure in the far field of the jet, R., and hence j? (), can be calculated 
directly. Assume that the sound wave is incident on the structure at an angle 8, as 
shown in Fig. 1. The velocity of the plane waves is c, the velocity of sound, and 
the velocity of the trace waves along the panel is c/sin@. The trace wavelength is 
assumed equal to A. 


Consider the case of the panel vibrating as a rigid body, giving y (a)=y (d)=1. 
When the trace wavelength A is very long compared with the span b there will be 
effectively a uniform pressure loading along the panel, in which case R., approaches 
unity. Therefore j?()=1 for b/A=0. It can readily be seen that when b/A=1 
pressures on one half of the beam will be balanced by suctions on the other half 
and j?(w)=0. The complete curve for the variation of j? (o) with b/A is given by 


A... 
7? (4 sin 


and is shown in Fig. 2. 


Consider now the case of the panel vibrating in a sinusoidal mode having n 
half-waves of wavelength b. In the special case when n=1 there will be a node at 
each support and all parts of the panel will have deflections in phase. The maximum 
excitation will occur when b/A=0, as in the previous case. When there are several 
semi-wavelengths in the panel, maximum excitation will occur when the pressure is 
everywhere in phase with the deflection, which is achieved when b/A=1. The 
variation of j? () with b/A is shown in Fig. 3 and is given by 


It can be seen from this that a typical structure, vibrating in a mode of n half-waves, 
acts as a wavelength filter. The highest response is obtained when the trace wave- 
length is equal to the mode wavelength. It can also be seen that the selectivity 
increases as the value of m increases. 


110 The Aeronautical Quarterly 


Ma} 


IS 


iS 


2 
OINT ACCEPTANCE) 


OINT ACCEPTANC 


EFFECT OF NOISE ON STRUCTURES 


10 


° 


° 


4 


WAVELENGTH RATIO b/A 
Ficure 2. Joint acceptance. Rigid panel. Ref. 2. 
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FiGureE 3. Joint acceptance. Flexible panel. Ref. 2. 
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2.4. EFFECTS OF THE SURROUNDING MEDIUM 


As the structure vibrates there will be forces exerted upon it by the surrounding 
air. In general these forces will be made up of two parts; one part in phase with 
the acceleration, and a second part in phase with the velocity. The component in 
phase with the acceleration causes a slight increase in the effective mass of the 
structure, but it is not considered that this will be appreciable in practical cases, 
The component in phase with velocity gives rise to a damping force due to the fact 
that the vibrating structure is radiating sound energy. This acoustic damping is 
considered in detail in Section 4.3. 


A third effect which has received little attention, but which may be important 
in special cases, is due to the effect of the air enclosed within the structure. 
Foxwell and Franklin®” have shown that at certain frequencies standing waves 
can be set up inside a cylinder, thus inhibiting movement of the cylinder walls. 


2.5. SIMPLIFICATIONS TO THE GENERAL THEORY 


It is shown in Section 4.3 that the damping of the normal modes is small. Thus 
the response will be primarily resonant, with only those forces having a frequency 
close to the natural frequency contributing appreciably to the excitation. As the 
power spectrum of the noise pressures is relatively flat, the power spectral density 
in the region of the natural frequency can be assumed to be constant. If it is 
assumed that j? (@) is also constant in this narrow band of frequencies, then the 
response equation can be simplified as follows :— 


WL (w,) Woo (w,) (w,) 


0 


0 


A? Wp, (0) 7? (for resonant response) 


= A? (,) 7? 
where 4, is the damping ratio. 


Therefore the mean square stress 


A? 
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As it has been assumed that the damping is small there will be little coupling 
between the modes of vibration. Thus it can be assumed that the mean square of 
the total response is equal to the sum of the mean squares of the responses in each 
separate mode. The mean square of the total stress is therefore 


s=n 


> Woo 7? (,) h,? 


2 3 
5,m,70, 


s=1 


Thus the r.m.s. stress level in each mode varies directly as the noise pressure level, 
and the joint acceptance, and varies inversely as the square root of the 
damping ratio. 


Having obtained the r.m.s. stress level and the loading spectrum, fatigue life 
estimates can be attempted. The effect of randomly varying stresses on fatigue 
life is not fully understood, but a first approximation can be made by using Miner’s 
“cumulative damage hypothesis.” Unfortunately this can be considerably in error 
for some materials. 


3. Noise Pressure Field 
3.1. INTRODUCTION 


It has been shown that both the power spectrum and the spatial correlation of 
the noise pressures are required to define the exciting forces acting on the structure. 
The mechanism of noise production by a jet efflux is so complicated that much of 
this information must be obtained experimentally. Lighthill has developed a 
theoretical analysis of jet noise in which the sources are assumed to be acoustic 
quadrupoles formed by the eddies in the turbulent mixing region. The qualitative 
results of this analysis are supported by tests on model and full-scale jets, but 
because the exact nature of the turbulence in the jet efflux is not yet known it is not 
possible to estimate quantitatively the amplitude and correlation of the associated 
noise pressures. Refs. 6 to 23 discuss the properties of the noise pressure field. 


3.2. POWER SPECTRUM OF NOISE PRESSURES 
3.2.1. Subsonic Jets: Near Field 


The main sources of noise are the eddies in the mixing region of the jet efflux. 
The typical size of these eddies is small at positions close to the exit nozzle but 
increases with increase in distance downstream. The smaller eddies produce noise 
in the higher frequency range and the larger eddies produce noise in the lower 
frequency range. The proportion of energy in the eddies which is converted into 
sound, however, is very small. The eddies will in general have two associated 
pressure components; one out of phase with velocity, which is the purely hydro- 
dynamic pressure proportional to 4V?, and a much smaller component, in phase 
with the velocity, which is radiating sound. Structure close to the jet is thus in- 
fluenced partly by hydrodynamic pressures and partly by sound pressures. A 
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further complication is due to the fact that noise components of different frequencies 
are emanating from different positions in the jet stream. Clearly this position is so 
complex that it is not possible to estimate analytically the pressure fluctuations on 
any structure in the near field. Measurements show that for a given jet the pressure 
amplitudes vary with jet velocity raised to a power of 2 close to the jet and a power 
of 4 at distances greater than about ten diameters from the jet centre line. 


3.2.2. Subsonic Jets: Far Field 


At some considerable distance away from the jet, the mixing region can be 
considered to form a point source and therefore the sound waves will be approxi- 
mately plane. The noise reaching a point in this far field thus takes the form of a 
randomly varying pressure wave. As would be expected from the composition of 
the noise source, this random wave has a power spectrum which is relatively flat 
over a wide range of frequencies. The noise field also shows marked directionality, 
the maximum total noise intensity occurring along a line making an angle of about 
30° to the jet centre line. 


Lighthill“* *” has developed the general theory of the production of noise by 
turbulence in the mixing region of a jet efflux. He shows that the noise is produced 
by a field of acoustic quadrupoles, whose strength depends on the degree of 
turbulence. Unfortunately the degree of turbulence is not known with any accuracy, 
but nevertheless the theory establishes the important parameters. It is found that 
the sound intensity varies typically as the jet velocity raised to the power eight, the 
so-called V* law; and as the square of the exit nozzle diameter. The sound pressure 
thus varies as V*. Experimental work on both model and full scale jets‘*-!*: '* 7°) 
shows that the general character of the noise field agrees with Lighthill’s 
predictions. 


3.2.3. Choked Jets 


The mechanism of noise production by a choked jet is usually complicated by 
the presence of shock waves in the efflux. Although a convergent-divergent nozzle 
can be designed to give shock-free flow at one particular operating condition, at 
other operating conditions shock waves will be formed. The interaction of eddies 
and hot spots with these shock waves gives off additional sound waves which 
increase the total noise output of the jet. Richards®® and Powell" have discussed 
these effects in detail. Vickers Supermarine@” have confirmed that the noise output 
of a full-scale jet is considerably increased when the jet is choked. For structural 
response estimations it will again be necessary to measure experimentally the power 
spectrum and the spatial correlation of the noise pressures. 


3.3. CORRELATION OF NOISE PRESSURES 
3.3.1. Space Correlation 


The next important property of the noise field from the structural response 
point of view has been shown to be the spatial correlation of the noise pressures 


114 The Aeronautical Quarterly 


CORRELATION COEFFICIENT Ry 


ovel 
at d 
forc 
pit, 
coel 


For 
the 
as § 
frec 
dep 
inci 
for 
ext 
for 
sid 
to | 


Co 
in 
2t 
reg 
no; 
dia 
ma 
the 
tio 


May 


= 


EFFECT OF NOISE ON STRUCTURES 
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CORRELATION COEFFICIENT Ry 


over the structure. This is a measure of the phase relationship between pressures 
at different points and hence of the effectiveness of the pressure field as an exciting 
force. Now as the pressure is fluctuating randomly the instantaneous pressures 
p(t,a) must be replaced by their time-averaged values. Thus the correlation 
coefficient in the narrow frequency band of central frequency w becomes 


Jim E he (t, a) p(t,a) at | 


0 


For discrete frequencies the pressure varies sinusoidally and hence the variation of 
the correlation coefficient with spatial separation takes the form of a cosine curve 
as shown in Fig. 4(a). In the near field of a jet efflux, however, there are no discrete 
frequencies and the noise field variation with distance is such that the inter- 
dependence between pressures at two points decreases as the spatial separation 
increases. Consequently the correlogram for the near field of a typical jet takes the 
form shown in Fig. 4(b) with R. falling off rapidly to zero. The full curve when 
extended to negative values of spatial separation is approximately symmetrical, and 
for practical purposes the distance between the first zero crossing points on either 
side of the datum represents the distance over which the pressures can be considered 
to be acting in phase. 


Full scale correlation measurements have been made by Callaghan, Howes and 
Coles®* on an axial flow turbine of 10,000 Ib. thrust rating. Analyses were made 
in eight frequency bands covering the range 100 c./s. to 1,000 c./s. in a region within 
2 to 26 nozzle diameters from the jet exit nozzle. These showed that the size of the 
region of positive correlation of the overall pressures varied from 0-385 to 0-93 
nozzle diameters. Work is being done at Southampton University on a 2 in. 
diameter model jet to study the mechanism underlying correlation patterns and to 
make comparisons with full-scale results. Franklin and Foxwell®* *® have shown 
that, for points close to and parallel to the boundary of the jet, the overall correla- 
tion can be deduced from the assumption that a fixed noise pattern is convected 
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downstream at a velocity equal to half the velocity in the centre of the jet. This is 
confirmed for full-scale jets but only holds for points close to the jet boundary. 


3.3.2. Auto-Correlation 


The auto-correlation coefficient is obtained by multiplying the signal by itself 
after a time delay 7 has been inserted into one channel. It is defined as 


lim E [ at| 
p (ty 
Although this is not required directly for structural response calculations it has the 
important property of being the Fourier transform of the power spectrum of the 
noise pressure. Thus the power spectral density of the noise pressure can be 
obtained from the relationship 


co 


W, cos w7 dr. 


0 


The variation of R- with time delay 7 can be measured and the Pegasus digital 
computer has been programmed to carry out the transform. This provides effec- 
tively a narrow band analysis and is thus a more accurate way of obtaining the 
power spectrum of the noise pressures than the more usual one-third octave filters. 


3.3.3. Measurement of Correlation Coefficients 


The measurement of correlation coefficients can be carried out most con- 
veniently by means of an electronic analogue computer. A general purpose 
correlator of this type, having a frequency range from 0 to 20 kc./s. has been con- 
structed at Southampton University. The input signals are amplified and pass to a 
quarter-squares multiplier. A bootstrap integrator integrates the product, which is 
then finally converted into digital form and registered on a Dekatron counter. The 
integration time can be any multiple of 10 sec. within the range 10 sec.-10 min. 
Integration times longer than 10 sec. are achieved by operating two integrators in 
sequence for 10 sec. periods. The 10 sec. standstill period of each integrator is 
used to transfer the result to the digital register and to correct for drift by means of 
a servo-mechanism. 


A time delay unit is now in operation with this correlator to enable auto- and 
cross-correlation measurements to be made. Using a frequency modulated system, 
the two signals are recorded on widely separated parallel tracks on a rotating 
magnetic drum. Corresponding reading heads are fixed some little distance round 
the circumference and one of these heads can be moved relative to the other to 
introduce a time delay in one of the channels. The two signals then pass to the 
multiplier and integrator in the normal way. 
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4. Mechanical Impedance of the Structure 


4.1. INTRODUCTION 


The determination of the mechanical impedance of the structure requires a 
knowledge of the normal modes and natural frequencies of vibration and also the 
damping associated with the normal modes. The analysis of the modes and fre- 
quencies of vibration of a cylindrical shell of uniform skin thickness has been 
presented by several authors. The general case of a stiffened shell of the type 
used in aircraft construction has not been fully analysed, although it is possible to 
make approximations to the lower modes of vibration. The inherent damping in 
an aircraft structure is usually low, but this can be increased by the addition of 
damping materials. As has been shown, any increase in damping reduces the 
amplitudes of vibration and hence the stress levels in the structure. 


4.2. NORMAL MODES OF VIBRATION 


42.1. Uniform Cylinder 


The most comprehensive analysis of the normal modes and natural frequencies 
of vibration of a cylindrical shell of uniform thickness is due to Arnold and 
Warburton®” *. Using the strain energy method and Lagrange’s equations, they 
analysed the flexural vibration characteristics of the shell and verified the results by 
tests on four steel cylinders of different dimensions. For the simply-supported 
cylinder the mode of vibration is assumed to be made up of sine waves longi- 
tudinally and circumferentially. The deflection of any surface element is written 


Axial displacement: u=A cos n¢ cos (mzx/1) 
Circumferential displacement: v=B sin n¢ sin (mrx/l) 
Radial displacement: w=C cos n¢ sin (mzx/1), 


where n is the number of circumferential waves and m is the number of longi- 
tudinal half-waves. The lowest root of the frequency equation gives the frequency 
of the mode having predominantly radial displacement. Curves are given by Arnold 
and Waiburton from which the variation in natural frequency with both n and m 
can be obtained for four different skin thicknesses. For a relatively thin shell it is 
found that, for a given value of m, the frequency decreases as n increases, which is 
contrary to what would have been expected. For thicker shells the reverse trend 
takes place. These results have been verified by tests on four steel cylinders having 
internal and external diameters of 3-748 in. and 3-950 in., and lengths ranging from 
384 in. to 15-63 in. 


4.2.2. Cylinder with Closely Spaced Stiffeners 

Miller®® has extended the work of Arnold and Warburton to apply to the case 
of a cylinder stiffened by frame rings and longerons. The spacing of the stiffening 
material is assumed to be small compared with the mode semi-wavelength in the 
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Ficure. 5. Vibration modes of a stiffened cylinder. Ref. 26. 


@ 


respective direction. It is then justifiable to assume that the mode shape takes the 
form of sine waves, as in the previous case. The strain energy and kinetic energy 
of vibration of the stiffening material are then added on to the energy of the uniform 
shell. This should give reasonably accurate results for the lower frequencies of a 
typical aircraft structure. Fig. 5 shows Miller’s estimate of the frequency curves 
for a stiffened cylinder of radius 60 in., skin thickness 0-048 in. and 15 ft. long. The 
stringer pitch is 8 in. and the frame pitch 20 in. In this case there are only 10 frame 
rings and 48 longitudinal stiffeners; thus the curves should give reasonable approxi- 
mations when 7 is less than 4 and n less than 10. These limiting values, however, 
need further investigation. 


4.2.3. Cylinder with Widely Spaced Stiffeners 


In cases where the stiffener spacing becomes comparable with the half wave- 
length in either circumferential or longitudinal directions it is not yet possible to 
calculate the natural frequencies. For the uniform cylinder, coupling terms 
between the modes were neglected. It seems probable, however, that these coupling 
terms assume considerable importance when the stiffener spacing is of the same 
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order as the mode semi-wavelengths. It is considered possible that the so-called 
panel and stringer modes, where one panel or stringer on the shell moves a greater 
amount than corresponding ones, can arise from these coupling effects. Some 
experimental evidence supports the idea of panel and stringer modes, but no 
thorough examination has been made. 


In practical aircraft structures, using stiffened cylinders, it is probable that 
over part of the range of frequencies excited by jet noise the modes can be deter- 
mined sufficiently accurately by assuming closely spaced stiffeners. Over the 
remainder of the range this approximation breaks down and the modes and 
frequencies May require experimental investigation. 


43. DAMPING OF THE NORMAL MODES 


4.3.1. Inherent Structural Damping 


Structural damping results from internal hysteresis of the material of the 
structure and from energy dissipation at the joints. The predominant effect is 
usually the damping due to the joints, as the energy loss due to hysteresis in the 
metal is very small. The damping due to energy dissipation at the joints depends 
on the normal pressure between the plates, due to forming the rivets, and also on 
the tangential force loading the joint. On loading the joint a point is reached when 
the frictional force between the plates is overcome and further increase in load 
causes slipping of the plates relative to one another. Most of the energy dissipation 
is due to this slipping. Thus a lightly loaded joint will not dissipate much energy, 
whereas a more heavily loaded joint is capable of appreciable energy dissipation. 


4.3.2. Methods of Increasing Structural Damping 


The amount of energy dissipated by joints can be increased by the use of a 
thin insert of damping material if an increase in the joint flexibility can be tolerated. 
This insert is subjected to a high shear strain, as the joint is loaded. Where the 
joint is made by adhesives, in the place of rivets, an improvement in the damping 
can be achieved by choosing an adhesive which is also a good energy absorber. 


With an existing structure the damping can be increased by the addition of 
damping materials, which are essentially layers possessing high energy absorption 
characteristics. Damping tapes have been used successfully by Vickers Super- 
marine®”, Alternatively “ anti-drumming ” compounds, developed originally for 
the motor industry, can be sprayed onto the inner surface of the structure. A high 
polymer type of plastic material is more efficient as an energy absorber than the 
older bitumen based materials and can be “blended” to give optimum damping 
over any particular range of temperatures and frequencies. This material is most 
effective when used on parts of the structure which carry the highest direct strains. 
Mead“? *») has investigated the damping properties of this type of material and 
Studied the effect of damping compounds in general on vibrations excited by 
jet noise. 
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Ficure 6. Variation in acoustic damping with frequency. 


4.3.3. Acoustic Damping 


Acoustic damping arises from the fact that as the structure vibrates it radiates 
sound and thus dissipates energy. The magnitude of the acoustic damping depends 
on the mode and frequency of vibration of the structure. Junger®® has analysed 
the case of a long cylinder radiating sound; and from this work a good approxima- 
tion to the acoustic damping can be made. The cylinder is considered to be vibrating 
with flexural waves both circumferentially and longitudinally and it is assumed that 
the cylinder is long compared with the longitudinal wavelength, so that three- 
dimensional radiation effects at the ends can be neglected. This may not be a 
severe restriction in practical cases of jet noise excitation. 


Typical values of the acoustic damping ratio have been computed for the 
cylinder discussed in Section 4.2.2. For modes having frequencies in the range 
200 c./s.-500 c./s. the damping ratio varied from 0-02 to 0:06. If it were possible 
to change the frequency of a given mode of vibration, the variation in damping ratio 
with frequency would take the form shown in Fig. 6. The curve has a marked peak 
and there is a “cut off” frequency below which no sound is radiated. Radia- 
tion does not take place when the structural longitudinal wavelength is less than 
the sound wavelength. In general it would not be possible to change the frequency 
of a given mode appreciably on an existing aircraft structure. This work does 
indicate, however, that in the design stage of a new aircraft it may be possible to 
achieve high acoustic damping ratios by a suitable choice of structure. 


5. Panel Tests 
5.1. INTRODUCTION 


In the early days of experience of fatigue failure attributable to the action of 
jet noise, many simple ad hoc tests were made in an effort to investigate the 
problem. Many of these were misleading because they failed to reproduce all the 
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TIME 
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(b) PANEL STRESS 


FIGURE 7. Resonant response of a panel to a random pressure wave. 


characteristics of the phenomenon. One extensive series of tests, being made 
primarily in the United States, involves the subjection of single panels to high 
intensity noise. The panels, having dimensions typical of aircraft fuselage panels, 
are usually mounted rigidly along their edges. The method of generalised harmonic 
analysis has been applied to the analysis of panel vibrations by Miles“. Subse- 
quently some agreement between theoretical and experimental results has been 
achieved, but the validity of extrapolating from these tests to complete aircraft 
structures has not been fully confirmed. 


5.2. THEORETICAL ANALYSIS 


Miles analyses the case of a panel vibrating in its fundamental mode due to 
excitation by random noise pressures. The wave of random noise is assumed to 
impinge normally on the panel; thus the correlation coefficient is equal to unity 
over the whole panel surface. The damping of this panel mode is small and, there- 
fore, the response is assumed to be essentially resonant. Thus the panel can be 
considered to respond only to those pressures having frequencies in the range 
©) +8w,, where w, is the panel natural frequency, and 4 is the damping ratio. The 
analysis is now simplified and the response takes the form of a vibration of 
frequency », but having an envelope which varies randomly, as shown in Fig. 7. 
The mean square value of the envelope amplitude is given by 


= (Wo) Sy” 


46 
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W, (,) is the power spectral density of the noise pressure at frequency «, and S, is 
the stress at the particular point due to uniform static pressure. 


This estimation of r.m.s. stress levels for panel tests gives good agreement 
with experimental results. The problem of fatigue failure under random loads, 
however, requires much more investigation before reliable estimates of fatigue lives 
can be made. 


5.3. EXPERIMENTAL RESULTS 


Lassiter, Hess and Hubbard®® have carried out tests on many panels to 
investigate the fatigue properties of different materials and different forms of con- 
struction. Flat panels 11 in.x13 in., varying in thickness from 0-032 in. to 
0-064 in., and curved panels having a radius of curvature of 48 in. were mounted 
rigidly onto one face of a supporting chamber. High intensity acoustic loading of 
the panels was obtained by means of a siren or a loudspeaker playing back noise 
recorded close to an actual jet efflux. The damping of the panels was measured 
by the resonant amplification method. Stress levels measured on the surface of the 
panel were found to be in good agreement with predictions based on Miles’s method 
of analysis. The tests were carried on until fatigue failure occurred, but no 
theoretical estimation of the fatigue life was achieved. 


Tests have been made by Vickers Supermarine on panels placed behind an 
actual jet engine, to investigate the fatigue properties of different types of construc- 
tion. It was found that the honeycomb sandwich type of construction is much 
superior to single skin panels. At de Hiavilland’s a limited ad hoc test on a 
complete aircraft®” indicated that a simplified analysis based on Miles’s work and 
assuming so-called panel and stringer modes could give reasonable agreement with 
measured r.m.s. stresses. More investigation is required before the general validity 
of this approach can be established. 


5.4. LIMITATIONS 


The main difficulty with the theoretical work of Miles and the experimental 
results of Lassiter, Hess and Hubbard and others is that it cannot be reliably used 
for the prediction of the behaviour of a complete aircraft structure until the validity 
of such an extrapolation has been established. In the first place the vibration 
characteristics of a stiffened cylindrical fuselage cannot be represented by a single 
panel. Secondly, the correlation pattern in the noise field cannot be reproduced by 
a loudspeaker or siren. A smaller discrepancy may arise from the fact that non- 
linear vibrations occur in panel tests, as a result of the membrane effect of direct 
stresses in the middle plane of the panel. Aircraft panels cannot be considered to 
be rigidly fixed along their edges; consequently non-linearity should not be present 
to the same extent. It can be concluded, therefore, that panel tests have value 
in checking the relative strength of different materials and forms of construction 
from a fatigue point of view; but they can give little indication of the magnitude 
of that fatigue life on an actual aircraft. 
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6. Fatigue Life Estimates 


6.1. GENERAL 


It has been seen that, having made certain reasonable assumptions, it is possible 
to estimate the r.m.s. stress levels in an aircraft structure but it still remains to 
provide some estimate of the fatigue life of the parts concerned. The effect of 
random loading on fatigue life has not been fully investigated and the fatigue 
properties of materials in the low stress range have received little attention in the 
past. One possible approach is to assume a normal distribution of stress and to 
use Miner’s “cumulative damage hypothesis ” to estimate the life under random 
loading conditions. Tests on simple specimens have indicated that this approach 
can be in error for low stress levels. Thus it should be emphasised that reliable 
fatigue life estimates cannot be made at the present stage of knowledge on the 
subject. In special cases it may be possible to make estimates based on full-scale 
test experience on a similar structure by using the theory to allow for the effect of 
structural changes. Much improvement can also be achieved by good detail design 
of attachments and joints. 


6.2. EFFECT OF RANDOM LOADING 


Freudenthal and Heller®® have carried out rotating-bending tests on plain 
polished 2024 and 7075 aluminium alloy specimens using random loading. Their 
results indicate that the fatigue life under random loading is lower than that pre- 
dicted by the cumulative damage hypothesis. The specimens were loaded by means 
of a lever system and an electro-magnet which produced, in a random order, six 
different stress levels. The stress levels were of an order associated with gust load 
conditions in an aircraft and not the low stress levels associated with jet noise 
vibrations. With gust loading, the load cycle is not‘completely random in sequence, 
as the highest loads are associated with other high loads in high level turbulence. 
It is not yet known whether or not this is true in the case of jet noise stresses. Thus 
the tests using a completely random sequence of loads may be severe. s 


The difference between the fatigue life achieved, and that predicted by the 
cumulative damage rule, increased as the r.m.s. stress level decreased for 2024 
aluminium alloy. At a characteristic random fatigue life of 10’ cycles the actual 
life had been reduced to only 10 per cent of the predicted value. For the 7075 
aluminium alloy there was more scatter of results, but the trend was towards a 
constant ratio of 0:7 between actual and predicted lives. These tests show that the 
effect of random loading is more severe than is predicted by the cumulative damage 
rule and also that there is considerable difference in the sensitivity of different 
materials to random loads. 


6.3. ESTIMATES OF FATIGUE LIFE USING THE CUMULATIVE DAMAGE RULE 


In attempting to estimate the fatigue life of a component subjected to randomly 
varying loads, it would be convenient to be able to express the loading in terms of 
an effective number of applications of the r.m.s. stress. If cyclic stresses, of 
magnitude equal to the r.m.s. stress, alone were applied, fatigue failure would occur 
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7 
in a time T= — sec., 


fo 


where — of applications of r.m.s. stress level required to produce 
failure 


f.=frequency of r.m.s. stress in c./s. 


If now the loading magnitude is varied randomly about this r.m-s. level, the 
life will be reduced, because stresses greater than r.m.s. value will produce more 
damage than stresses lower than r.m.s. value. To allow for this Powell® introduces 
the factor 7 such that 


T= e. for random loading. 


e 


Using the cumulative damage rule, assuming a normal distribution of stresses, and 
a net peak representation of the loading due to Rice“, Powell shows that 


dy 


0 
where 


N=number of applications of stress o to produce 
failure 


1/2 


w?W, dw 
1 


and 


where w, (w) is the power spectral density of stress at frequency . 


Thus, when the fatigue curve of the material is known, 7 and f, can be estimated 
for the case of random loading. 


However, much more experimental investigation is required to justify this 
approach. 
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Note on the Forces that Act near the Centre 
and the Tips of a Swept-Back Wing 


D. G. HURLEY 


(Aeronautical Research Laboratories, Melbourne) 


SUMMARY: The momentum equation is used to calculate the drag that acts 
near the tip of a semi-infinite swept-forward wing of constant chord that is 
at zero incidence to a uniform stream of incompressible, inviscid fluid. The 
drag is given as a function of wing section and angle of sweep and is shown to 
be unaffected by shaping the tip. The results may be used to calculate the drag 
that acts near the centre and the thrusts that act near the tips of a swept-back 
wing, provided that its aspect ratio is sufficiently great for the regions where 
the different forces act to be separated. Some results due to Neumark and to 
Kiichemann and Weber suggest that this will be so provided that the aspect 
ratio is greater than about two. The results are combined with some due to these 
authors to estimate the spanwise extent of the forces. It is found that this 
extent decreases as the angle of sweep is increased and that it depends quite 
markedly on the wing section. 


1. Introduction 


The low speed potential flow about untapered non-lifting swept wings has 
been considered by Neumark” and by Kiichemann and Weber”. One of their 
findings was that wing sections near the plane of symmetry experience quite large 
drag forces, whereas sections near the tip experience balancing thrusts. They also 
found that when the aspect ratio is greater than about two the regions where these 
forces act are separated by ones where the flow is very nearly the same as that 
about the corresponding swept wing of infinite span—a flow in which the force per 
unit span is zero. 


These findings suggested that it should be possible to use either the energy or 
the momentum equation to calculate the drag and thrust forces in these cases and 
this paper shows how this can be done. 


The results are useful, since the approximate method used by the authors 
teferred to is more suited to calculating the drag that acts at the centre section 
rather than the total drag that acts near the centre section, which is given quite 
simply by the present method. 


It is convenient to introduce the method by considering the flow about a semi- 
infinite swept-forward wing. 
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NOTATION 


a@ unit vector in direction of V, 
c chord measured normal to leading edge 
C contour of wing section described in positive sense 


Cp=D/(4pV 7c”), drag coefficient of tip of semi-infinite swept-forward 
wing cut off by plane normal to leading edge 


Cy =D’/(pV?c”), drag coefficient of tip of semi-infinite swept-forward 
wing cut off by vertical plane containing V, 


C, = d./(4pV ,*c), drag coefficient of wing section at centre of swept wing 


D drag of tip of semi-infinite swept-forward wing cut off by plane 
normal to leading edge 


D’ drag of tip of semi-infinite swept-forward wing cut off by vertical 
plane containing V, 


D, drag of triangular portion of wing (see Fig. 3 on p. 136) 


d, drag force per unit span that acts on wing section at centre of 
swept wing 


F resultant force on wing tip 

Lone half edge length of cube (see Fig. 1) 

m source strength 

p fluid pressure 

P. free stream pressure 

S surfaces of cube excluding S, (see Fig. 1) 

S, cross-sectional area of wing section (see Fig. 1) 
t thickness of wing 


Uz,U, perturbation velocity components of two-dimensional motion past 
wing section when speed at infinity is unity 


V.V_ fluid velocity and speed, respectively 
V.,V,. free stream velocity and speed, respectively 
V.,V,,V. components of fluid velocity 
Uz, Vy,Uz perturbation velocity components (see equation (1)) 


rectangular 


co-ordinates leading edge from tip, 


z vertically upwards (see Fig. 1) 


x normal to leading edge, 
Xx, 


/ 


x’ co-ordinate measured in direction of V, 
B angle of yaw (see Fig. 6) 
A angle of sweep 
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w doublet strength (see equation (6)) 
p fluid density 
¢, velocity potential: u,= —09¢,/0x, — 09,/0z 


vy, stream function corresponding to ¢, 


2, The Semi-Infinite Swept-Forward Wing 


Consider the low speed potential flow about a semi-infinite untapered wing of 
symmetrical section which is swept forward through an angle \ and which is at zero 
incidence to a stream of speed V, (see Fig. 1). 


Let Oxyz be a set of rectangular axes with O at the tip and Oy along the 
leading edge. Ox lies in the chordal plane of the wing and Oz is vertically 
upwards. Let V., V,, V. be the components of the fluid velocities and let v,, vy, vs 
be the perturbation velocity components, so that 


V,=V, cosrA+v,z 
V,=2, 


The momentum equation is applied to the fluid which lies within the large 
cubic control surface S shown in Fig. 1, and this gives 


F=-({pas—| . Q 
KY 


Ss 


{SOMETRIC VIEW 
FiGurE 1. Semi-infinite swept-forward wing, showing control surface and notation. 
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Here F is the force acting on the portion of the wing inside the cube, d$ is a vector 
element of surface area in the direction of the outward normal to S and V is the 
fluid velocity. 


Now, by Bernoulli’s equation, 


=p. sind+v,)?+v7}, from equation (1) 
= — pV, (vz cos A+, sin A)—4p (v.27 +. . ‘ eve 


Substitution into equation (2), and taking the scalar product with a unit vector 
a in the direction of V,, gives 


pa. | | ds+p (v, cos A+, sin A) Via. dS+ 
s 
+4p (v2? | V.a(V.4d9), 


where D is the component of F in the direction of V,. 
Since V.a=V,cosA+V, sind 


=V,+v,cosA+v,sinA, from equation (1), 


D=-p,a. | (vzcosA+v, sind) V..dS+ 
Ss Ss 


+4p | | (V,+v2 cos A+¥v, sin A) (V. dS) 
Ss 


=—p,a. | dS+p | cos A+v, sin A) (V,—V).dS + 
s s 


| | | ; 4 
s 


Now since the surface S excludes the surface S, where the wing protrudes 
through the cube, 


a. | sina, 


so that the first term on the right hand side of equation (4) is p,S, sinA. The last 
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term on the right hand side is zero because f V .dS is the total flux of fluid 


across S, which is zero. 


Thus equation (4) becomes 


D=p,S, SinA+p (v, cos A-+v, sin A)(V,—V).dS+4p +v,?+,7)a.d8. 

(5) 

The limit of this equation as L —> © is considered, where L is half the edge 

length of the cube. It is assumed that the velocity perturbations at large distances 

from the wing can be considered as arising from a semi-infinite line doublet of 


strength » lying along the positive y-axis and a source of strength m at the origin. 
With these assumptions it can be shown that 


“t (x? (x? +y?+27)j/2 (x? + 27) (x? + y? (x7 + y? + 
my 
(x? + y? + 27)3/2 (x*+y? + 27)3/2 

(x? + 27)? (x? + y? + /2 (x? + 27) (x? + y? /? (x?+y? + 


(6) 


This equation shows that v,, v, and v, are O(1/L*) on all the faces of the cube 
other than the face y= +L. On this face v, is O(1/L*) but v, and v, are not. 
Also, on this face, ; 


(V,-—V).dS = —v, dS, 


so that the limit as L —> 00 of equation (5) is 


D=p,S, sind +4p | (v,? +0,2) a. dS 


oo 


= p,S, sind +4 sin A j @2+0)dS, . . . 


co 


where the region of integration extends over the whole of the xz-plane, at y= +0. 
It is evident* that, as y—> +00, v, and v, approach the values for the two- 
dimensional flow about the wing section when the speed at infinity is V,cosA. 
‘The replacement of the wing by a semi-infinite line doublet (and a source), as has been done, 
is nearly equivalent to replacing it by a semi-infinite circular cylinder. On the face y= +L 


the values of v, and v, given by equation (6) differ from the two-dimensional values for a 
circular cylinder only by quantities which are O(1/L?). 
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If u, and u, denote these two-dimensional velocity perturbations when the Speed at 
infinity is unity, then, as y—> +00, 


Uz—> u,V, Cosa, 


so that equation (7) becomes 


D=p,S, sind +4pV,? sin A cos? A | — 
Now +u,’) dS= § . . 
co Cc 


(Ref. 3, p. 66), where ¢, and ¥, are the velocity potential and stream function 
corresponding to the perturbation velocity components u, and u, and C is the 
contour of the wing section. Thus equation (8) becomes 


_ D__ pS,sind 
* J ody, . . (10) 


This equation is the same as equation (9) of Ref. 4, which was derived using 
the energy equation. 


3. Physical Interpretation of Equation (8) 


If the whole of the wing surface were subject to a uniform pressure p,, then 
the component of the force acting on the wing in the direction of V, would be 
pS, sin A, which is the first term on the right hand side of equation (8). It is thus 
seen that this term represents a buoyancy type force and is of little interest. 


To obtain a physical interpretation of the second term on the right hand side, 
it is convenient to multiply the equation by V, and to consider the motion as arising 
from the wing’s being advanced through fluid which is at rest at large distances from 
it. The left hand side is then equal to the rate at which work must be done to 
produce the motion. 


A cubic control surface is again introduced, but it is now considered to be at 
rest, so that the length of wing inside the cube is increasing at a rate V, sind. The 
second term on the right hand side of equation (8) (multiplied by V,) can be written 


{V, sin A} x { 4 [ f [uM cos A)? +(u.V, cos ry | dS } ; 


oo 


The second factor in this expression is the kinetic energy of the fluid per unit span 
of wing where it protrudes from the cube (the y velocity component being now 
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FIGURE 2. Variation of tip drag of swept-forward wing with angle of sweep. 


zero). Thus the expression is the rate at which the fluid inside the fixed cube is 
gaining kinetic energy and, as remarked previously, the left hand side of the 
equation represents the rate at which work must be done to generate the motion. 


4, Dependence of C, on Angle of Sweep and Wing Section 


Equation (10) shows that, apart from a buoyancy type force, Cp varies as 
sind cos? A for all wing sections and this function is shown in Fig. 2. It is zero for 
A=0° and 90° and has a maximum when A=35° 16’. 


The effect of differing wing sections was investigated by evaluating 


j 


for 2 number of standard aerofoil sections, using the ordinates and surface velocities 
given in Ref. 5. 


¢, was obtained by numerical integration from 


s 


Where s is the distance measured around the surface of the aerofoil from its 
leading edge. 
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is given by ¥,=-z(x), . ; 


where z=z (x) is the aerofoil profile (Ref. 3, p. 83). 
1 


was evaluated numerically, using the values of ¢, and ¥, given by equations (11) 
and (12). The results are given in Table I, which also contains the exact value for 
an elliptic section, (~/4)/(t/c)? (Ref. 3, p. 85). 


Approximate calculations were also made for the biconvex section 


z(x)=2— 0<x<c, i «u » 


which is considered extensively in Refs. 1 and 2. 


Equation (2.4) of Ref. 2 is, with c=1 and V,=1, 


1 


0 
which gives u {2+(x-1) log (+*)} (15) 


To the same degree of approximation as is involved in equation (14), ¢, may be 
evaluated from 


which gives, using equation (15), 
{ x-+xlog, x—(1—-) log, (1 x)- x? log. x + (1 x)? log, (1 - x) } 


It then follows from equations (12), (13) and (16) that 
1 sry’ 
Cc 


a value which is included in Table I. 
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TABLE I 
1 
Wing Section a dv, c J 
NACA. 63-010 0-00667 066) 
65-010 000673 
66-010 0:00693 0-693 
65-006 000243 0-676 
65-008 0:00429 0-670 
Ili ( * =0°785 
Ellipse =) 57 
i 2 
Biconvex (equation (13)) * _09-637 
™\e 


The first four aerofoils listed in Table I are each 10 per cent thick, but differ in 
the location of the point of maximum thickness. The table shows that this has little 
eflect on the drag. The third, fifth and sixth aerofoils have similar thickness 
distributions but different values of t/c. The third column shows that the drag 
varies almost exactly as (t/c)*, which is also seen to be the case for elliptic and 
biconvex sections. The table also shows that, for a given value of t/c, an ellipse 
gives a larger value of the drag than the other sections. 


5. Wing Tips Cut Off in Direction of V, Instead of Normal to 
Leading Edge 


To enable comparisons to be made with the results of Refs. 1 and 2, it is 
necessary to consider the case when the wing tip is cut off by a vertical plane which 
contains V, and not by a plane which is normal to the leading edge, as has been 
the case so far. 


If D’ is the drag that acts in the former case, then 
D’=D+D,, . ‘ . (18) 


where D, is the drag that acts on the triangular portion of wing shown in Fig. 3. 


From this figure, 


c 
dz 
D,=2 cos 2 px tana dx 


0 


=2 sind px dx. , (19) 
0 
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FIGURE 3. 


When the triangular portion of wing under consideration is sufficiently remote 
from the tip for the surface pressures to be the same as for a swept wing of infinite 
aspect ratio, the use of Bernoulli’s equation leads to 


D.=2sind | { 3 u2)} x 
0 
q 
0 


since dx=S,, 
x 


the cross-sectional area of the profile, so that, from equations (8) and (18), 


c 
0 


(20) 


It is shown in the Appendix that the linear theory approximation to the second 
term on the right hand side of this equation is zero for all profiles. It can also bk 
shown*, using the results of Ref. 6, p. 160, for example, that the term is exactly ze10 
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for the case of elliptic profiles. The term was also evaluated numerically for the 
N.A.C.A. 65-010 profile, using the data of Ref. 5, and once again the result was, to 
within the accuracy of the calculations, zero. 


Although no proof that it is exactly zero in the general case was found, the 
preceding results suggest that it can be neglected in practice, so that equation (20) 
can be written in coefficient form 


__D’___ sind cos? A j 
2 


using equation (9). 


6. Application of Results to Swept Wings of Finite Aspect Ratio and 
Other Cases 


Consider a swept wing of large aspect ratio, as shown in Fig. 4. Now if the 
aspect ratio is sufficiently large, then the flow near the tips will not be influenced by 
the distortion to the flow produced by the kinked centre section, and the preceding 
analysis may be used (with the sign of A changed) to calculate the net force in the 
direction of V, that acts on each of the tips. Equation (21) shows that each force 
will be a thrust and, since the total force that acts on the complete wing must be 
zero, it follows that there will be a drag force acting on the central portion which, 
in coefficient form, will be of magnitude 2Cy., with Cp. given by equation (21). 


An idea of how large the aspect ratio must be for this separation of the forces 
to apply can be obtained from some calculations carried out by Neumark”. He 


FIGURE 4. 
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FicurE 5, Swept wing differing from one of constant chord section and sweep only near centre, 


considered a biconvex wing of aspect ratio 2 having 53° 8’ sweep and found that, 
even for this small aspect ratio, the velocity field inboard of mid-semi-span was 
nearly the same as when the aspect ratio is infinite. It thus appears that if the 
aspect ratio is greater than about 2 it is legitimate to treat the central and tip 
portions of the wing separately and to apply the preceding analysis. Some further 
comments on this are made in Section 8. 


The preceding results may also be applied to a wing of infinite aspect ratio 
whose shape departs from one of constant chord, sweep and section only near the 
centre (see Fig. 5). For the disturbances to the “infinite aspect ratio flow” at 
sections of the wing remote from the centre may be represented by a certain 
distribution of doublets, and these disturbances will die out sufficiently rapidly for 
the preceding analysis to apply. It is thus seen that the total drag force that acts 
near the centre section of a swept wing of large aspect ratio cannot be reduced by 
changing the sections near the centre. This remark, of course, applies only to low 
speed flow and not to transonic flow. 


The present method can also be used to treat the case of a wing which is swept 
through an angle A and also yawed through an angle f (Fig. 6). The thrusts acting 
on the port and starboard tips will be given by equation (21) with A replaced by 
—A+f and -A-—f respectively. The drag force acting near the centre will be 
numerically equal to the sum of the two tip thrusts. 


9. The Results of Neumark” and of Kiichemann and Weber”? 


These authors have calculated the flow that occurs near the centre of a swept 
wing of infinite aspect ratio by a method in which the wing is replaced by a system 
of kinked sources. Kiichemann and Weber obtained the following expression fot 
the drag coefficient of the central section (see equations (3.13) and (3.26) of Ref. 2). 
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+B 
FIGURE 6. 
+ sin az x’ 


is the wing profile in the direction of V,. 


To obtain the value of 2Cp, the total drag that acts near the centre, it is 
necessary to know the values of the local drag at sections of the wing other than the 
central one and then to integrate across the span. This involves lengthy numerical 
calculations which appear to have been carried out only for the biconvex section 
of equation (13). Figure 12 of Ref. 2 gives results when A=45° and Fig. 19 of 
Ref. 1 gives results when A=53° 8’. 


For the biconvex section, equation (22) gives 


) cos? A log. x)+0 o(t ) ‘ (23) 


and by numerical integration of the two figures just referred to it is found that 


2Cyp = 0-42 (t/c)? for A=45° 
and =0-°37 (t/c)? for A=53° 8%’. 
The values of 2Cp, given by the present method (equations (21) and (17)) are 
0-45 (t/c)? for A=45° 


and 0:37(t/c? for A=53° 8, 


which are seen to agree reasonably well with the results of the other authors. 
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2D'= TOTAL DRAG ACTING NEAR CENTRE OF SWEPT WING 
d, =DRAG PER UNIT SPAN ACTING AT CENTRE SECTION 


0-20 
0-15 
D’ 
x 
0-10 


30 60 90 
FiGurE 7. Estimates of the spanwise extent of the central drag force for wing having biconvex 
sections. 


8. The Spanwise Extent of the Drag Force that Acts Near the Centre 
of a Swept Wing 


As remarked in the previous section, it is only for swept wings of biconvex 
section that the method of Refs. 1 and 2 has been used to calculate the flow at 
sections other than the centre one, so that knowledge of the spanwise extent of the 
region where the drag acts is very meagre. Additional information can be obtained 
by calculating the ratio D’/d,., which has the dimensions of length and is a measure 
of the spanwise extent of the region where the drag acts. 


For the biconvex section, equations (17), (21) and (23) give 


D3 _csinA 
d. 1+sinA 
Be 1—sinA 


and the main term is seen to be independent of t/c. 


Figure 7 gives the values of D’/(cd.) and it is seen that this decreases as 4 is 
increased. 


For wings having elliptic sections, equation (21) and Table I give 


2 
D’=4pV sin A cos? A (+) 
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2D'= TOTAL DRAG FORCE ACTING NEAR CENTRE OF SWEPT WING 
a, == DRAG PER UNIT SPAN ACTING AT CENTRE SECTION 


0°20 
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D’ 
0-05 
30 60 90 


FicurE 8. Estimates of the spanwise extent of the central drag force for wings 
having elliptic sections. 


and equation (3.28) of Ref. 2 gives* 


=f 2 2 
d.=4pV, A log. 


where k? = 1 —[(t cos A)/c]* and K (k) and E (k) are the complete elliptic integrals of 
the first and second kind, respectively. 


Thus, for elliptic sections, 


D_ c sinA k? 
d, 16, /1+sind\ K(kK)-E(K)’ 
log. 
1—sinA 


Figure 8 gives the values of D’/(cd.) for various values of t/c. It is seen that 
D’/(cd.) increases as t/c is increased, unlike the case when the wing section is 
biconvex, when it was independent of t/c. It is also seen that D’/(cd.) decreases 
as \ is increased in a similar way to when the wing has a biconvex section. 


*It is noted that as r/c > 0 this equation gives the surprising result that d,/(t/c)? >, so that 
the spanwise extent of the drag force would appear to tend to zero. 
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Appendix 


c 
EVALUATION OF | Qu+ut-+ug) x (dz/dx) dx FOR ARBITRARY PROFILES 
0 
UsING LINEAR THEORY 


Introducing the approximations of linear theory and taking c=1 the integral becomes 


dz 
=2 

I Zax 

0 
dz _ 
1 

and 


where P denotes that the principal value of the integral is taken. 
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FORCES ON A SWEPT-BACK WING 


1 1 

0 0 


1 


°l—x 


0 
and, from equation (24), 


xo (x) [o ()—o (x)] 2 
1=2 {| dtdx+ ~ | xo (x) log. dx 
00 0 
0 0 
11 
1,= = (x) dt dx (26) 
1 
and /,= 2 | x0" (x) log, 
0 


If do /dx exists, I, is proper, so that the order of integration may be changed to give 


1 


00 


and interchanging the symbols x and t gives 


2 
2 to (t) 
|| 

00 


Thus, from equation (26), 


1 
2 
1,+1,= 4 + dx—= | at ds 
0 0 00 


(t) dt dx, 
0 


0 0 0 


since the value of [xo? (x)— to? is unchanged by interchanging x and 1. 


May 1959 143 


ly, 

Lift 
igs 
on, 
an, 
ion 
ty, 

| 

1eS 


D. HURLEY 


Thus, using equation (24), 


xo? 


| 
— dt dx+ — ~ | P| 


0 0 0 0 


and is thus zero, since z=0 at x=0 and 1. 


Using a theorem due to Hardy", the order of integration in the third term may be 
changed, provided that do /dx is continuous, so that it becomes 


1 1 1 
0 0 


0 0 


to (t) log, 
0 


The second term on the right hand side of equation (27) is also equal to this 
expression, so that equation (27) becomes 


1 
2 1— 
I,+1,= (t) log, dt 
0 


=-—I,, from equation (26). 


Thus, from equation (25), 7=0. 
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On the Buckling of Certain Optimum Plate 
Structures with Linearly Varying Thickness 


(27) 


E. H. MANSFIELD, M.A., Sc.D., A.F.R.Ae.S. 
(Royal Aircraft Establishment) 


SUMMARY: This paper is concerned with the buckling under uniform 
longitudinal compression of a variety of structures composed of plates whose 
thickness tapers linearly to zero across the section. Such structures include the 
angle of Fig. 1, the strut of cruciform section of Fig. 2 and the simply-supported 
strip of Fig. 3. For given cross-sectional area and overall dimensions (e.g. 
length of arm) the sections with linearly varying thickness achieve a greater 
buckling load (assuming that local buckling, rather than Euler buckling, is the 
criterion) than sections with any other smooth variation of thickness. These 
particular sections are therefore optimum sections and, even if they may not 
be used in practice, provide a convenient yardstick for purposes of comparison. 
The buckling loads are considerably greater than those for the corresponding 
this “constant thickness” sections. 

NOTATION 


Ox,Oy Cartesian axes, Ox measured lengthwise, Oy measured across the width 
from the simply-supported edge 


a_ length of plate 

b width of singly tapered plate 

t, maximum thickness of tapered plate (at y=0) 
t thickness of plate 

E,G Young’s modulus and shear modulus 

vy Poisson’s ratio 

D flexural rigidity, Er’ / {12 (1 —’)} 

a=t,/b 

oc, uniform compressive buckling stress 


w displacement normal to plane of plate 
n integer 


A, arbitrary constant. 


Received February 1959 


erly May 1959 145 


y be 


E. H. MANSFIELD 


NS 


FIGURE 1. FIGURE 2. 


1. Method of Solution 


The basic problem is the determination of the buckling stress of the tapered 
plate shown in Fig. 4 subject to conditions of simple-support along the ends x=0, a 
and along the edge y=0, and zero moment and shear along the edge y=b. Because 
of the variable flexural rigidity of the plate the differential equation") contains 
derivatives of D: — 


dD 0 | 0?w 
4 2 
DV w)+ dy? Yaa + dy? +o,t =(). (1) 


Et,? y\3 


@D Et (t,\? 


and equation (1) can therefore be written in the form 


A glance at equation (3) shows that if 
the last term vanishes, and if 


the first three terms vanish. There are an infinite number of solutions of equation 
(5) satisfying all the boundary conditions and the general solution can be expressed 
as follows in the form of an infinite series with arbitrary coefficients. 
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BUCKLING OF PLATE STRUCTURES 


FIGURE 3. FIGURE 4. 


w= = A, sin (nzx/a) sinh (nzy/a). (6) 


Equations (4) and (6) represent the complete solution, which might also have 
been deduced from the particular approach of Ref. 2, based on energy 
considerations. It is to be noted that, because of the vanishing of D and dD/dy 
along the free edge, the conditions there are automatically satisfied by any smoothly 
varying deflection function. 


The fact that the buckling stress is proportional to the shear modulus G, rather 
than to E/(1—v’), is at first surprising, until it is seen that in any of the buckling 
modes the sum of the principal curvatures is zero, so that conditions of pure twist 
exist everywhere. The fact that there is an infinite number of possible buckling 
modes occurring at the same value of the compressive stress implies that the 
structure is an optimum one; any small variation of thickness designed to raise the 
buckling stress in certain modes will lower the buckling stress in others. 


If the edge of the plate is clamped, or elastically restrained against rotation, the 
plate with a linearly varying thickness is no longer the optimum one; it will, in fact, 
buckle at the same stress as for the simply-supported case—the buckles being 
initially confined to a vanishingly narrow strip along the free edge : — 


w =~ Lim A, sin (ntx/a) sinh (nty/a). (7) 


n->0o 


2. Comparison with Constant Thickness Sections 


For purposes of comparison with the angle of Fig. 1 and the cruciform strut of 
Fig. 2 it is sufficient to consider a single plate simply-supported along one edge and 
free along the other. The thickness of the constant thickness plate of equal weight 
is 4¢, and the buckling stress is therefore given by 


for the tapered plate. 
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It follows from equations (8) and (9) for moderate ratios of a: b that the tapered 
plate, and hence the angle and the cruciform strut in Figs. 1 and 2, will buckle at a 
stress nearly four times greater than one of constant thickness. However, for very 
short plates (a < 0-95) the constant thickness plate is the better; this is because the 
influence of the supported edge is diminished and the plate tends to buckle as an 
“Euler plate”. Similarly, for very long struts, in which Euler buckling of the 
complete strut is possible, the constant thickness strut may again be better. 


If a> 2b, the simply-supported constant thickness plate, of thickness 41, and 
width 26, buckles at a stress given by 


It follows from equations (9) and (10) that the doubly tapered plate will buckle 
at a stress 24(1—v)/x? times that of the constant thickness plate, and it will be 
noted that 


24 (1—v) 


or=1-82 if v=0-25 


These ratios may be compared with the maximum value of 1-40 obtained by 
Capey™) for a stepped cross section with a thin central strip. The doubly tapered 
plate is, of course, ruled out on practical considerations because of the vanishing 
thickness at the centre, resulting in zero shear rigidity in the plane of the plate. 


FIGURE 5. 


However, a localised increase in thickness over this region, such as indicated in 
Fig. 5, could give the necessary shear rigidity and strength without seriously altering 
the buckling stress for the “lower” modes, the buckling stress for the “higher” 
modes being considerably increased. 
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Flow in a Viscous Trailing Vortex 


B. G. NEWMAN, M.A., Ph.D., A.F.R.Ae.S. 


(Department of Engineering, University of Cambridge) 


SUMMARY: The equations of motion for an isolated laminar viscous vortex 
at moderate to large Reynolds numbers are linearised, by assuming that both 
the rotational velocity and the deficit of longitudinal velocity are small 
compared with that in the free stream. The rotational motion and the 
longitudinal motion may then be superimposed and solutions are readily 
obtained for each. If the vortex is generated by a body with profile drag it 
is predicted that the deficit of longitudinal velocity will be positive, which is 
in agreement with experimental observation. Further details of the solution 
and its relation to the flow in real vortices are discussed; and the theory is 
compared with some measurements in a turbulent vortex. 


1. Introduction 


When fluid flows past a lifting wing trailing vorticity is produced which, at some 
distance downstream, eventually concentrates into two trailing vortices. At large 
distances from the centre of each vortex the rotational or swirl velocity associated 
with each is very closely equal to that of a free irrotational vortex. This fact is 
exploited when determining the effect of vortices at large distances as, for example, 
in the determination of ground effect and wind-tunnel interference. Near the centre 
of each vortex, however, the distribution of rotational velocity progressively deviates 
from that of a free vortex and at the centre the velocity is zero instead of infinite. 
A better approximation to the trailing vortex is then provided by the Rankine 
combined vortex consisting of a core of fluid rotating like a solid body, surrounded 
by a free vortex. In this model of the flow the radius of the forced vortex core 
should increase as the flow proceeds downstream and the vorticity diffuses outwards. 
The rotational velocity at any particular radius within the core decreases until, very 
far downstream, the velocity is negligible and the vortex has completely decayed. 
However, an inviscid model of this kind does not enable the diffusion of vorticity 
and decay of the vortex to be predicted. Furthermore it is implicitly assumed that 
the longitudinal or axial velocity within the vortex is constant. This does not appear 
to be true in general; for example, Hilton” has measured a reduced longitudinal 
velocity near the centre of a trailing vortex. 
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The diffusion of vorticity and decay of a vortex are important in some practical 
applications. Andrews and Kraft have both shown that the flow disturbance 
several span lengths downstream of an aircraft is mainly associated with the 
rotational flow in the trailing vortices and thus a knowledge of their decay is of 
interest when assessing the behaviour of a second pursuing aircraft. Another 
example is afforded by the use of small, low aspect ratio fins as vortex generators 
to prevent boundary-layer separation from a wing or other surface. The trailing 
vortices re-energise the boundary layer and the downstream extent of their effective. 
ness depends on the extent to which the rotational velocity within the vortices 
is maintained. 


In this paper a small perturbation analysis of a single trailing vortex is 
presented. The vortex is formed in a viscous incompressible fluid and is assumed 
to be laminar. The equations of motion and continuity are linearised by assuming 
that both the rotational velocity and the deficit of longitudinal velocity within the 
vortex are small compared with the undisturbed stream velocity. As a result of the 
linearisation the equations for the rotational and longitudinal motion separate, and 
simple solutions are obtained for each. A similar type of solution has been given 
by Gortler® for the flow in a jet with weak swirl, although it differs from the 
present analysis in that the variation of static pressure associated with the swirl 
is neglected. 


NOTATION 
r,6,z cylindrical polar co-ordinates: Oz is along the axis of the vortex 

u__ radial velocity component 
v rotational velocity component 

w’ longitudinal velocity component parallel to Oz 

W longitudinal velocity in the free stream at large r 
w = W —w’, longitudinal velocity defect 
static pressure 

p,, Static pressure in the free stream 
total pressure 

P,, total pressure in the free stream 
p density of the fluid 
v kinematic viscosity of the fluid 
n eddy viscosity when the flow is turbulent 

I’ circulation round the vortex at large r 

D_ drag of body producing the vortex 

D, profile drag, value of D as I‘—> 0 

D, induced drag 

r,=2z(v/Wz)'/?, radius of the equivalent forced vortex 


150 The Aeronautical Quarterly 


radial 


rotati 


axial 


and 


undi: 


Choo 
and 
tudin 

whe 

Ole 

May 


FLOW IN A VISCOUS TRAILING VORTEX 


v,=(1'/4xz) (Wz/¥)'/”, rotational velocity at the edge of the forced vortex 
R=r/r,, non-dimensional radius 
V =v/v,, non-dimensional rotational velocity 


@ measured upwash angle within the vortex in a plane parallel to the 
trailing edge and the free stream velocity, ~ tan-' (v/W). 


2. Analysis of a Laminar Viscous Trailing Vortex 


The flow is assumed to be steady, axi-symmetric, incompressible and laminar. 
Choosing cylindrical polar co-ordinates r, 6 and z, where Oz is the axis of symmetry, 
and denoting the radial velocity by u, the rotational velocity by v and the longi- 
tudinal velocity by w’, the equations of motion are 


du. v1 2 [ 
Ow’ Ow’ 1 
ial: +w = -— +vV’? ‘ 
axial : +w > +vV?w, (3) 
0? 1 a 0? 
where V?= ort ag 
and the continuity equation is (ru)+ ==@), ‘ (4) 


Putting the axial velocity w’.=W-w, where W is the axial velocity in the 
undisturbed fluid at large r, it is assumed that 


v and w_ are small compared with W, i.e. O (<) if W is O(1), 
u_ is very small compared with W, i.e. O (¢’). 
Thus, if z is O (1), the continuity equation (4) indicates that r is O (¢). 


The radial equation (1) therefore becomes 


@) €&) © (v) (vy) (ve”)_ 


ou ou ou 1 Oop 1 ou, 4 


Where the order of magnitude of the various terms is indicated. 


If the Reynolds number of the main flow Wz/y is large and O (1/<”), i.e. v is 
0 (¢*), then the radial equation simplifies to equation (5) which follows : — 
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also indicating that (1/p) @p/dr is O 


The rotational equation (2) becomes 


ov ov uv 1 dv 


and thus simplifies to 


0 
The axial equation (3) becomes 
oz Oz or? r Or | 
and hence this equation simplifies to 
ow sis Ow 


It will be noticed that equations (6) and (7) are now linear and the variables v and 
w have been separated. 


The continuity equation (4) becomes 


1 a ow 


By the nature of the approximations we are examining the vortex some distance 
downstream of its origin. Hence it is sufficient to assume that the vortex is suddenly 
generated at z=0 as a free vortex of circulation I’. Far downstream the vortex 
finally decays until all the perturbation velocities u, v and w are once again Zero. 
The boundary conditions to be satisfied by equations (5) to (8) are thus 


z = 0, v=I'/(2tr) and w=0 except at the singular point r =0, 
z > 0, v and w— 0 when r is large (though still small compared with 2), 


zZ—> oo, v and w—> 0 for all r. 
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FLOW IN A VISCOUS TRAILING VORTEX 


The solution of equation (6) which satisfies these boundary conditions is 
identical with that given by Lamb“ for the development with time of a two- 
dimensional viscous vortex if the time ¢ is replaced by z/W, and it has previously 
been given in this form by Squire“ : — 


T v2) 

A —Wr'/(4v2) 
The solution of (7) is w= e ‘ ‘ (10) 


where A is a constant, as yet undetermined. 


The distribution of radial velocity u is given by equation (8):— 


Ar — Wr? /(4vz) 


973 © (11) 


uaz 

The constant A may be determined by considerations of momentum. The 

rotational and the axial-radial flows are independent of one another and thus 

equations (10) and (11) represent the distributions of w and wu for all values of the 

vortex strength I’. Thus if D is the drag of the body which is producing the vortex 
it is apparent that the limit of D as '—> 0*, D, say, is such that 


| p (W —w) w2tr dr= pWw 2trdr, 
0 
since w is small compared with W. 


D, 


Hence A= 


The static pressure distribution is dependent only on the rotational velocity v 
ind is given by equations (5) and (9). 


For large r, v= a “free” vortex 


’ 


8rvz 


and for small r, v= , a “forced” vortex. 


These “free” and “forced” vortices meet at 
y 1/2 T Wz 1/2 
(pz) and - 


‘In proceeding to this limit it is assumed that the wake flow remains axi-symmetric. 
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r, may be thought of as the radius of the vortex core of an equivalent Rankine 
vortex (see Fig. 1). 


It is convenient to rewrite equation (9) with v and r made non-dimensional in crag 
terms of v, and r,. sia 
Thus V= l-e where V=v/v, and R=r/r,. 
The maximum value of V is 0-639 and occurs at R=1-12. 
The distribution of V is shown in Fig. 1. 
The variation of static pressure is given by equation (5):— 
1 op 1 | 
—, +e-7"’]. T 
pv,” OR Rl | in the 
Hence, integrating from R to ©, 
Po-P_ 1 _ i(—R? ] i(—2R? ] 
in? +Hi(-R)|+]  +2Bi(-2R%) |, 
Hower 
smalln 
where Ei (— x)= -{ —— dt, one of the exponential integrals, and p,, is the static | contrit 
signific 
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pressure in the free stream. The limiting value of (p,,— p)/(4pv,") when R —> 0 is 
log. 2, and the variation of (p,,— p)/(4pv,”) with R is shown in Fig. 2. 


The distribution of total pressure P through the vortex is dependent on both the 
| totational and the longitudinal motion. 


where P.,, is the total pressure in the free stream. 


Thus, if w and v are of the same order of smallness, the deficit of total pressure 
|in the vortex is associated mainly with the longitudinal motion and 


D, e v2) 
2nxpyWz 


However, it is not essential to the analysis that v and w are of the same order of 
smallness, but merely that both are small compared with W. Thus the second-order 
static | contribution to the total pressure variation from the rotational motion may be 
‘ignificant when the longitudinal velocity deficit is very small. 


arterly 
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2 
This is given by = on 


=2 [Ei 


which is also shown in Fig. 2. 


3. Discussion 


The present solution is based on the assumption that the rotational velocity 
and the defect of longitudinal velocity within the vortex are both small compared 
with the free stream velocity. It is therefore clear from equations (9) and (10) that 
the solution is invalid in the immediate vicinity of the origin O where the vortex is 
generated, but that it becomes increasingly accurate as the flow proceeds 
downstream. 


Certain details of the solution are worth noting. The radius of the equivalent 
forced vortex core r,=2z(Wz/v)~*/? and is proportional to z'/*. As expected from 
the nature of the approximations, the growth of the vortex core is similar to that of 
a laminar boundary layer in zero pressure gradient. The maximum rotational 
velocity Vmax =0°639 (I'/4xz) (Wz/¥)'/? occurs at a radius 1-12 r, and at this radius 
the static pressure rise from the centre of the vortex is just over half the total rise 
attained at large radii.* 


The maximum deficit of total pressure occurs at the centre of the vortex and is 
given by 


D, 

If the vortex is assumed to be one of those generated by a wing with an elliptic 
spanwise distribution of circulation, then the vortex drag D,=47pI™ and equation 
(12) may be written 


} 


the two terms being respectively the contributions associated with the longitudinal 
and the rotational motion. Since D, is usually of the same order as, or less than, 
D,, the contribution from the longitudinal motion is usually greater than that from 
the rotational motion if, as seems probable, most of the wake leaving the wing is 
gathered into the two trailing vortices, so that 2D, represents a large proportion of 
the total profile drag of the wing.t This tendency is noted in the measurements of 
Fage and Simmons in the wake behind a wing. 


*Compared with the Rankine vortex for which exactly half the pressure rise occurs at the 
edge of the forced vortex core. 


{Note that the total drag of the wing is expressed as 2D, + Dj. 
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FLOW IN A VISCOUS TRAILING VORTEX 


Since the effective profile drag is usually positive, there will usually be a deficit 
of longitudinal velocity in the trailing vortices generated by a wing. Such a deficit 
has been measured by both Hilton® and Templin and it has also recently been 
found by Mowforth within the vortex generated by a half wing mounted on the wall 
of a wind tunnel (full details of these measurements are given in the Appendix). This 
conclusion will presumably apply generally to vortices generated by a closed body. 
However it is easy to imagine a situation in which the deficit might be negative, for 
example behind a jet-flapped wing. 


When an attempt is made to relate the present solution in greater detail with 
the flow in a real trailing vortex, two questions arise. Can an isolated trailing vortex 
of infinite extent exist and will the flow within it be effectively laminar? 


The deficit of longitudinal momentum within the vortex has been equated to the 
profile drag of the vortex-producing body on the assumption that all the wake is 
gathered into the vortex. The remaining induced drag D; must clearly balance the 
net pressure force on a control volume surrounding the body. 


co 
Thus D,= | (p.,— p) 2tr dr. 
0 


At larger, p,,—p Oc + , and hence D, —> 


In other words, the drag of a body producing an isolated trailing vortex of 
infinite extent would be infinite and therefore such a vortex cannot strictly exist. A 
finite wing, of course, produces two trailing vortices, and the rotational velocities 
associated with each cancel at large radii so that D;, is finite. A single trailing vortex 
of finite extent is readily produced in a duct“, in which case the flow field does not 
extend to infinity and the drag is once again finite. As usual the effect of the duct 
boundaries may be determined by the introduction of a suitable series of image 
vortices. Within the limits of the present small perturbation analysis, the equation 
for the rotational velocity, and hence the trailing vorticity, is linear and is thus 
unaltered by the superposition of an irrotational flow field. Hence the velocities 
associated with the various real and imaginary vortices may be superimposed 
vectorially for as long as the effectively rotational parts of the vortices (r << 2-5r,) do 
not overlap.* Hence, with this restriction, the present solution may be applied 
directly by superposition for the trailing vortices produced by a wing or for a single 
trailing vortex surrounded by boundaries. 


It is more difficult to reach a satisfactory conclusion concerning the second 
question. The flow which is gathered into the trailing vortex from the wake will in 
general be turbulent and not laminar. However, the circulation increases with 
radius and thus there is a tendency for the turbulent. flow to be attenuated and 


*It is interesting to note that this conclusion is implicit in Prandtl’s determination of the radius 
of the forced vortex core within the two Rankine vortices trailing behind a wing.) 
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possibly completely damped out, particularly near the core of the vortex. The high 
degree of stability for a rotational flow in which the circulation increases with radius 
was first pointed out by Rayleigh and was noted by Taylor“” for the flow between 
rotating cylinders. It was also found that when the speed of rotation was sufficient 
for the flow to become unstable and turbulent the associated turbulent stress was 
nevertheless small. It is therefore expected that the eddy viscosity associated with a 
turbulent vortex will be small compared with that usually associated with free 
turbulent flow. 


Following Gortler® and Squire’, it is attractive to apply the present laminar 
solution to a turbulent vortex by simply replacing the kinematic viscosity by an eddy 
viscosity which is constant at each downstream station*. The present measurements, 
which were made some three semi-spans downstream of a rectangular half-wing, are 
therefore compared with the theory in this way. The results are shown in Figs. 3 
and 4, where the measured values of downwash angle, Pitot and static pressures are 
compared directly with the theory for various values of the eddy viscosity. Full 
details of the measurements and the method of reducing them will be found in the 
Appendix. The results indicate that an eddy viscosity of approximately 10v is 
appropriate for the deficit of dynamic pressure associated with the longitudinal 
motion near the centre of the vortex, with somewhat higher values farther out, 
although it should be noted that the vortex is no longer axi-symmetric in this region. 


*This is considered to be slightly more convenient than the usual procedure of adding an eddy 
viscosity to the kinematic viscosity. 
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For the rotational motion the independent measurements of static pressure and 
downwash angle are in agreement in giving the following values for the eddy 
viscosity near the centre of the vortex : — 


12 in. station, 7=50v 
16 in. station, 7=30v 
and farther out »=100v approximately in both cases. 


It should be noted that these results are qualitatively in agreement with 
Rayleigh’s ideas in that the eddy viscosity is least in those regions where | dI"/ dr | 
is greatest. It is clear therefore that the flow cannot be strictly described in terms 
of a constant eddy viscosity at each downstream station. However, for the central 
part of the vortex 7 is approximately constant and for some engineering purposes 
the deviation from this curve at larger radii may be unimportant. Perhaps it should 
be noted that the deviation is comparable with that found for similar analyses of 
turbulent jets. The idea that a constant eddy viscosity gives a useful engineering 
approximation to the real fluid flow will therefore be pursued. 


The reduction of 7 in the downstream direction for the rotational motion is 
probably associated with the fact that the vortex was not fully rolled up. Templin® 
has noted that when this is so the apparent viscosity will be large. The vorticity at 
any radius is greater than that in the fully rolled-up state: thus the vorticity appears 
to have diffused to larger radii and this leads to a fictitiously high value of the 
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eddy viscosity. The fall in computed Cp, with downstream distance (see the 
Appendix) and the lack of complete axial symmetry support the hypothesis that the 
vortex was not fully rolled up. 


From dimensional considerations it would be expected that n/v would be a 
function of I'/v and Wz/¥, probably becoming independent of Wz/v when it is large, 
Mowforth’s values for 7/v, and also those of Templin and those of Titchener and 
Taylor-Russell*, are presented in Fig. 5, where the values of Wz/v are shown in 
brackets. 


T 
® TEMPLIN 


MOWFORTH 
on © TITCHENER @ TAYLOR-RUSSELL 
i 
69x10") (FIGURES IN PARENTHESES ARE THE 
1 «0 ROTATIONAL VALUES OF Wz/v, WHERE z IS MEASURED 


(93n10°) FROM THE CHORD OF THE VORTEX 
GENERATOR) 
20 oy | 


x10° 
i73*65x10°) | 
LONGITUDIN L (6-4+6 9x1) @7-2-101x! *) 
5 


10 
Vv 


Ficure 5. 


These few results seem to be consistent in indicating that /v falls to a constant value 
with increasing Wz/v and, for large values of Wz/v, tends to unity as I'/v exceeds 
about 8 x 10*. However this curve should only be accepted tentatively. /v would 
hardly be expected to remain unity for very high values of I'/v. For instance, 
Timme* has noted that the flow in the cross-stream vortices forming the Karman 
vortex street behind a circular cylinder is effectively laminar if I‘/v is Jess than 
2-5 x 10° and that above this Reynolds number the vortices become turbulent, with 
an apparent eddy viscosity roughly ten times the kinematic viscosity. Furthermore 
free turbulent flows usually become independent of viscosity, which should give »/T 
constant at large Wz/» (Ref. 7). Clearly more measurements are needed before 
Fig. 5 can be accepted with any confidence. 


4. Conclusions 


A small perturbation analysis for a laminar viscous vortex has been given. This 
indicates that there is usually a deficit of longitudinal momentum within the vortex 
which is directly associated with the profile drag of the vortex generator, while the 
cross-stream variation of static pressure associated with the rotational motion 
directly balances the induced drag. The solution has been compared with measured 
turbulent vortices by replacing the kinematic viscosity by a constant eddy viscosity. 
The agreement is only fair, but it is noted that the theory may be useful for some 
engineering purposes. 


*Templin’s vortex was measured very near the trailing edge of a rectangular wing of aspect 
ratio 1/4; that of Titchener and Taylor-Russell was generated in the centre of a circular pipe. 
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FLOW IN A VISCOUS TRAILING VORTEX 


Appendix 
DETAILS OF E. MOWFORTH’S MEASUREMENTS IN A TURRULENT TRAILING VORTEX 


The vortex was produced by an untwisted rectangular wing, of length (equivalent 
semi-span) 4-20 in. and chord 1-56 in., which was mounted as a half-model in the middle 
of the 20 in. wall of the 74 in. X 20 in. low speed wind tunnel at Cambridge University. 
The wing, which was uncambered, was set at 6° incidence and had an aerofoil section 
of 184 per cent thickness/chord ratio with maximum thickness at 25 per cent chord. The 
wind speed was approximately 100 ft./sec. and was adjusted throughout the experiments 
to maintain a constant Reynolds number per foot at the wing of W/v=6-60x 10°. The 
value of W/v increased slightly down the tunnel, because of the growth of the boundary 
layer on the tunnel walls; the values of Wz/v are noted on Figs. 3 and 4. A trailing 
vortex was thus formed in the centre of the tunnel and was measured at two stations, 
z=12 in. and z=16 in., downstream of the quarter-chord point. The following measure- 
ments were made : — 


(i) Total pressure P, by means of a 0-028 in. diameter Pitot tube with lips sharpened 
internally. This tube measured the true total pressure for angles of yaw up to 
12° and had an error of only one per cent of the dynamic pressure at 20° yaw 
(note that the measured yaw did not exceed 11°). The tube was traversed 
across the vortex, both parallel and perpendicular to the trailing edge. 


(ii) Angle of upwash ¢, with a null-reading Conrad yawmeter made from two 
0-028 in. diameter tubes. The yawmeter was traversed through the centre of the 
vortex and parallel to the trailing edge. 


(iii) Static pressure p, with a 0-028 in. diameter static tube aligned with the local 
flow direction as indicated by the yawmeter. This was also traversed parallel 
to the trailing edge. 


The centre of the vortex was determined as the point of minimum total pressure 
from a preliminary survey with the Pitot tube. The measurements of downwash and 
static pressure have been corrected for small interference of the traversing gear. The 
downwash has also been corrected for the interference of the nearest tunnel walls 
(i.e. those 74 in. apart). This correction was small and did not exceed 0-1°, which was 
the accuracy of measurement. The associated correction to static pressure was negligible. 


As the vortex is not strictly axi-symmetric, the measurements are compared directly 
with the theory rather than in terms of a mean rotational velocity and a mean longitudinal 
velocity defect. The theoretical curves are presented for various values of the eddy 
viscosity 7. For the downwash and static pressure measurements the theory has been 
matched with the measured value of gr at large r. The values at both stations are in 
agreement in giving I'/v=1-97 X10‘. Using lifting-line theory for a rectangular wing at 
6° incidence, this gives a two-dimensional lift-curve slope of approximately 5-2, which is 
a plausible value, bearing in mind the low Reynolds number of the experiment and the 
fact that about 4 in. of the wing was immersed in tunnel-wall boundary layer at the wing 
root. The mean curve of measured static pressure has been subtracted from the Pitot 
pressure to give the results shown in the top curves of Fig. 3. The theory has been 
matched to the value of [(P,,—P)—(p,,— p)]/(ZpW?) at the centre of the vortex. Assuming 
that the eddy viscosity 1=10v, the measurements give 


Cy, =0-00183 at 12 in., Cp,=0-00190 at 16 in. 


These values are very low but it should be noted that much of the deficit of 
longitudinal momentum is found ‘1 the outer part of the vortex, when the eddy viscosity 
is apparently greater than 10v and where, in any case, the flow is not strictly 
axi-symmetric. It also is interesting to note the slight increase of C,, with distance 
downstream. Discussion of these measurements is given in the main body of the paper. 
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A Contribution to the Theory of Flow 
Stability in Multi-Stage Axial Compressors 


A. D. S. CARTER, B.Sc., A.M.I.Mech.E., A.F.R.Ae.S. 


(Royal Military College of Science, Shrivenham) 


SuMMARY: Previous theories on the surging of axial compressors have been 
based on a quasi-static theory, whereas a finite rate of growth of the circulation 
around each blade must occur in practice. An examination of this factor is 
made in the paper, studying the stability with respect to small changes in speed 
instead of mass flow as in the normal approach. Effects of capacity are 
neglected. It is concluded that an instability occurs in the last stage of the 
compressor, even when that stage is itself unstalled, which could give rise to 
surge. It has been found necessary to introduce one empirical factor into the 
analysis but, when this is given a value compatible with other experimental 
and theoretical evidence, reasonably good agreement with experimental results 
is obtained. 


1. Introduction 


Although the surging of multi-stage axial compressors is a feature of paramount 
importance in gas turbine design and operation, the problem has received compara- 
tively little experimental or theoretical attention. Experimental investigations are 
considerably hampered by safety considerations and usually confined to the simple 
determination of the surge line. Theoretical work is probably limited by the 
absence of basic test data on which to formulate any hypothesis, while the conflict 
of information available on surge behaviour may be a contributory factor. 


The basic mechanism of surging has been expounded by Kearton™ and Den 
Hartog’’. A more formal mathematical treatment of this basic concept has been 
given in Reference 3. It is well known, however, that the treatment does not accord 
with experience, although the resulting conclusion that surge will occur at peak 
pressure ratio is often used as a practical guide. A much more advanced study was 
made by Pearson and Bowmer"?. However, for the average compressor this gives 
a surge line which does not depart to any significant extent from that based on peak 
pressure ratio. 


A fundamental criticism of the previous theoretical work is that a unique 
stage or compressor performance characteristic is assumed. Stability is usually 
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studied by the method of small displacements in the operating parameters. Any 
lag in change of performance which would accompany such a small change in 
operating conditions should really be included, but has always so far been 
neglected (Pearson and Bowmer quote lag as one of the major factors which could 
modify the results obtained by their theory). Furthermore, although the stability 
of operation of a compressor with respect to mass flow changes is studied, stability 
with respect to speed changes appears to be neglected. These aspects of the problem 


are taken up in this paper. 


NOTATION 


annulus area 


constant in equation for temperature rise coefficient 


F 

k constant 
P pressure 


ratio of circulation during transient to circulation with constant axis‘ 


velocity 
R gas constant 
s blade pitch 
T temperature 
U_ blade speed 
Vs fluid velocity 
W mass flow 
a fluid angle measured from axial direction 
AP pressure rise 
AT temperature rise 
I’ circulation 
n stage efficiency 
Q work done factor 


0 at exit from previous blade row 
1,2 at entry to, and exit from, blade row 
I,II at entry to, and exit from, compressor 

a axial component 

t constant throttle conditions 

s stage value 

w_ whirl component 

7 value during transient speed increase 
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FicurE 1. Flow conditions in a simple single-stage system. 
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2. General Theory 
2.1. STABILITY OF A STAGE 

Consider in the first instance the stability of a stage which forms part of a 
multi-stage compressor. The conventional approach involves the arbitrary displace. 
ment of one of the operating parameters, say mass flow. If the system is to be 
stable, then it can easily be shown (see, for example, any of the references previously 


quoted) that 


where AP is the pressure rise across the stage, V, is the fluid axial velocity, suffix U 
denotes constant blade speed (or r.p.m.) and suffix t denotes constant throttle setting, 


Figure 1(a) shows some typical compressor stage characteristics and also a 
typical constant throttle line. It is seen that the flow is stable according to the 
criterion (1). It is also seen from Fig. 1(a) that this can be written in the less 
familiar form 


This condition can also be shown to be necessary by the physical nature of the 
problem. Suppose the speed increased from U, to U,, as defined in Fig. 1, then 
the compressor would be giving more pressure rise than is required by the throttle. 
To achieve this pressure rise the compressor speed would drop, by virtue of the 
excess load on the blades, thus making for stability. If the condition (2) were not 
Satisfied the pressure rise at speed U, would be less than that required by the 
throttle and so the speed would have to rise, making for instability, as shown in 
Fig. 1(b) for example. This condition locates a possible point of instability at the 
peak pressure point as usually assumed for practical purposes. It is shown in 
Appendix I that (1) and (2) are only equivalent if the stage is unstalled. In that 
case the efficiency and the fluid outlet angle can be taken as independent of the fluid 
axial velocity. It can then be shown (as in Appendix I) that the conditions (1) and (2) 
for stability can be written 


0a, 0a, 


Condition (3) can also be derived directly by assuming that a small change of mass 
flow takes place while the speed remains constant, and condition (4) can be derived 
by assuming that a small change in speed takes place while the mass flow remains 
constant (as given in detail in Appendix I). A unique stage performance in terms 
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of mass flow and speed has to be assumed in both cases. Now it does seem 
reasonable that any small adventitious mass flow changes will occur in such a 
manner as to satisfy the internal aerodynamics of the stage. However, speed 
changes can be quite arbitrary and illustrate a fallacy in the simple stability 
argument. Suppose the speed increases by a quantity 6U from U to (U + 6U). 


At speed U the circulation will be given by 


where |* is the circulation around the blade, and V,, and V., are respectively the 
whirl velocity at the inlet to, and the exit from, the blade row. 


It can then be shown (as in Appendix I) that 


al 


The increase of speed (8U), supposedly taking place while the axial velocity remains 
constant (due to inertia) in the stability argument, must therefore be accompanied by 
an increase of circulation round the blades given by 


An increase of circulation cannot occur any more than an increase in axial 
velocity. This is more clearly demonstrable if the speed increase 5U is supposed 
to take place suddenly or instantaneously, a postulate usually made in the orthodox 
stability argument. Likewise, both axial velocity and circulation cannot remain 
constant. This difficulty can only be resolved by supposing that the sudden or 
instantaneous speed increase takes place in finite time, and that the axial velocity 
and circulation both change in a mutually compatible manner in that time. To take 
the most general case it can be assumed that, following a rapid speed increase, the 
circulation changes by p times the change which would be obtained if the axial 
velocity could have remained constant, where p is some as yet unknown factor. 


Then _ = ps . (8) 


where the suffix 7 denotes the value during the transient speed increase 6U. 
From this relation it follows, as indicated in Appendix I, that 


{? tan z, —(1 — p) tan z, + tan a, \ cos? 2, (9) 


A similar expression for (62,/0V.), could be obtained mathematically by 
supposing a sudden mass flow increment, in which case the compressor speed would 
have to adjust itself to accommodate the rate of growth of circulation. 
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Ficure 2. Stability of a stage assuming a finite rate of growth of circulation. 


An examination of isolated aerofoil data and calculations (see Ref. 5 for 
example) suggests that the circulation would grow very rapidly at first, but would 
take some time to reach its ultimate value. A value for p somewhere between 0:5 
and 0:85 might be expected. Failure of the circulation to attain its full value for 
constant axial velocity displaces the momentary operating point from R to S, as 
shown in Fig. 2(a). It has been pointed out to the author that this régime is only 
possible physically if there is a capacity between the compressor and the throttle. 
Such an assumption is necessarily made although, in studying the stability, any other 
effects of capacity are neglected. For the example shown in Fig. 2(a) the system 
is still stable, but the displacement may be so large, or the throttle line so orientated, 
that instability results, as shown in Fig. 2(b). The condition for stability which 
replaces (4) is 


0a, 0a, 

au), for stability . . (10) 

tan z,—(1 —p) tana, +tan2, | cos’ a, 
(si tan a, + tan a, Vz 


If the axial velocity could remain constant during the transient speed increase 
(i.e. p=1), the condition for stability would become 


If the circulation remains constant during the transient speed increase (i.e. p=), 
the condition for stability is given by (13), which follows. 
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Oa, _ cos’ a, 3 


The latter is clearly the more rigorous condition. As nearly always, reality probably 
lies between the extremes, and the condition for stage stability will be applied to a 
multi-stage compressor retaining p as an unknown factor. 


a - 
% 
Ww 
w 
= 
a < 
a 
z 
SPEED U SPEED U 
a * 
% 
” 
« 4 
a 
MIDDLE 
z 
& STAGE = 
SPEED U SPEED U 
a 
<4 
< 
x 
a 
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FicurE 3. Operating conditions for various representative stages in a multi-stage compressor. 
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Ficure 4. Variation of stability derivative (0a,/U), through a multi-stage compressor. 


2.2. STAGE OPERATING CONDITIONS AND STABILITY OF A MULTI-STAGE MACHINE 


It will be appreciated that, for any stage in a multi-stage compressor, the constant 
throttle operating line will be dictated by the characteristics of the throttle itself plus 
the characteristics of the remaining stages of the compressor. The throttle 
and those stages could be regarded as a composite throttle, whose character- 
istic is different from the simple example, and different for each stage in the 
compressor. If any stage becomes unstable it is assumed that the compressor as a 
whole surges. The variation in operating conditions with speed has been given in 
some detail in Fig. 18 of Ref. 6. In Fig. 3 these have been replotted using the 
AP, U and z,, U co-ordinates. It can easily be seen from Fig. 3 that (0a,/0U); 
will increase from the first to the last stage, as illustrated in Fig. 4. Since the condi- 
tion for stability in the absence of stalling postulates that this must be less than a 
given quantity, it follows that critical conditions will always be reached initially in 
the last stage of the compressor, assuming that it is a normally matched machine 
with similar blading throughout. Surge will therefore originate in the last stage of 
the compressor if the type of instability discussed in this paper is the operative one 
in practice. The conditions pertaining to the last stage before surge will correspond 
to those illustrated in Fig. 2(a), the surge point being transition from the conditions 
in Fig. 2(a) to those in Fig. 2(d). 


It may be noted that the last stage has been singled out as the critical stage for 
a compressor in which no stage stalls. In an actual compressor the first stage of 
stages will stall, and the simple analysis adopted in this paper will no longer apply. 
However, Pearson and Bowmer") have shown that the orientation of the constant 
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throttle and constant axial velocity characteristics of these stages usually satisfies 
stability when a unique characteristic is employed. Assuming finite rate of growth 
of circulation (and separation) will make stalled stages more stable, so no generality 
appears to be lost. Attention is therefore concentrated on the last stage. 


It will be convenient to consider in the first instance pressure ratios greater than 
2:1. This implies a choked throttle orifice, so that if losses in the outlet ducting 
system are neglected the constant throttle lines can be defined analytically by 


W/ (Tx) 


Ir 


where W is the mass flow and T;; and P; are respectively the outlet temperature and 
pressure from the compressor. 


It is shown in Appendix II that in these circumstances 


By substituting this value of (02, /0U), in (11) the condition of stability becomes 
oAT 
(Gr ), tan + tan a, & 


It can be seen that if the simple approach is followed, based on an instantaneous 
change of circulation corresponding to a constant axial velocity during the transient 
speed change (i.e. p= 1-0), then the condition for stability is 


This is an obvious condition which could have been written down directly, but it 
supports the method of analysis adopted in this paper when it is applied to a simple 
case. Examination of any compressor test result shows that this condition is not 
nearly approached at surge, and stability will be investigated for values of p less 
than unity. 


Along any constant throttle line it is possible to write 


where F is a constant depending on the throttle setting, and f (U) is some function of 
U which represents the effect of speed on the matching. With this assumption 
(e@AT/0U), can be calculated. Substituting this value in (16), the condition for 
stability of the last stage, and hence of the compressor as a whole, becomes (see 
Appendix II) that given in (19), which follows. 
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Va (l-p) Tu 1 
stage {AT (U)} (tan + tan (19) 


where f’ (U) is the derivative of f (U) with respect to U. By means of this condition 
the stability of any point on a characteristic can be tested. The surge point is the 
highest flow at which instability is found. The expression would not apply to 
pressure ratios below 2:1, and even at this pressure ratio ducting losses may be 
playing a significant role. At medium speeds, when the throttle is just not choked, 
no simple analysis seems possible, but at very low speeds a constant throttle line 
can be represented by 


AP=kV> . ‘ (20) 


In this case the condition for stability becomes 


U last stage (tan Ao + tan (approximately). (21) 


For the range of values for z, and 2, used in practice, and if p is limited to the range 
of values already suggested (and more particularly if limited to the practical values 
determined later), then the flow coefficient for instability will be very small. This 
does not exclude any other form of flow breakdown, such as separation and rotating 
stall, which could occur at higher flows and in the practical range. In fact this does 
occur, so that according to this hypothesis there must be a fundamental difference in 
the mechanism of “‘surge”’ at high speeds compared with that at low speeds or with 
incompressible flow. 
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FiGurRE 5. Variation of non-dimensional temperature rise with speed at constant throttle setting. 
Throttle setting chosen as mean value at which surge was recorded in tests. 
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3. Comparison with Test Results and Discussion 


Very few compressor test results have been published in sufficient detail to 
facilitate a direct comparison of theory and experiment. One such set of results is 
given in Ref. 7 and another set is given (in less detail) on the same compressor, but 
with some blading modifications, in Ref. 8. The same value of p should therefore 
apply to both cases and enable a direct check to be made on the known factors in 
the condition for stability. 


In making comparisons it must be noted that the tests were carried out at wide 
speed increments (though not at all wide by normal standards), so that the calcula- 
tion of the function f (U), and more particularly its derivative f’ (U), is open to some 
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FiGuRE 7. Comparison of calculated and measured surge points on basis of flow coefficient to 
last stage for compressor of Ref. 7. 


error. The function f(U), together with the test points by which it is defined, has 
been plotted for both compressors in Fig. 5, and the derivatives are shown in Fig. 6 
to illustrate the errors involved. In Figs. 7 and 8 a direct comparison is made 
between calculated and measured flow coefficients at entry to the last stage 
corresponding to surge. 


At high speeds the calculated value of the flow coefficient lies in the stall region of 
the stage. Consequently (02, /0U) , will not be constant, as assumed in the analysis. 
Numerical work is impossible, although it can be seen qualitatively that the value of 
the flow coefficient at which instability is first reached will be increased by stalling. 
Consequently it is assumed in Figs. 7 and 8 that stalling of the last stage itself 
represents a surge limit, an assumption which is confirmed by many results on 
50 per cent reaction compressors. 


It will be seen that good agreement between the calculated and measured flow 
coefficients is obtained at other speeds. The general shape of the curves would lend 
confirmation to the theoretical criterion. Particular attention is drawn to the hump 
in both the curves in Fig. 8 between 8,000 and 8,500 r.p.m. One or two more test 
points would have confirmed the correlation, but the values at 8,500 r.p.m. are so 
exceptional that the chances of coincidence are small. In these circumstances it is 
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Figure 8. Comparison of calculated and measured surge points on basis of flow coefficient to 
last stage for compressor of Ref. 8. 


claimed that the test results confirm the hypothesis advanced in this paper. Closer 
agreement cannot be expected until more detailed test results are available on the 
unsurged performance (or more reliable methods of calculating it). The function 
f (U) could then be estimated accurately. In this connection it is interesting to note 
that agreement is closest where the function f (U) is most accurately defined. Such 
regions include the hump already referred to and the speed range 7,000 to 
8,000 r.p.m. of the comparison made in Fig. 7. 


Clearly the general level of the calculated curves is directly proportional to the 
empirical factor (1—p) and in that respect any correlation is meaningless. The 
value of 0-83 taken for p suggests that the circulation has in fact changed consider- 
ably, perhaps somewhat more than might at first be expected. However, in Fig. 9 
the rate of growth of circulation round an isolated aerofoil (taken from Ref. 5) has 
been plotted and the relevant value of p indicated. It shows a remarkable coinci- 
dence with the flattening of the curve, that is with the value of the circulation which 
would be maintained for some time during a transient. 


The theory presented in this paper differs from previous work in two respects. 
The major difference is the assumption of a finite rate of growth of circulation during 
the small speed increase by which the stability is studied. The second difference is 
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FIGURE 9. Range of possible values for p based on isolated aerofoil data (see Ref. 5, Fig. 181). 


one of means rather than principle and involves the examination of the effect of 
small speed changes rather than small mass flow changes. It is essentially a one- 
dimensional theory. 


A consequential difference is that, apart from the fact that changes in work 
done and efficiency, particularly in the earlier stages of a machine, may considerably 
alter the shape of the constant throttle lines, stalling does not play any direct réle in 
the phenomenon. The instability itself originates, at least at certain speeds, in the 
unstalled last stage, whereas the stalled first stages have hitherto been held 
responsible. At these speeds agreement between theory and experiment is best. 


It is most probable that a complete theory would have to include capacity 
inertia and so on, along the lines followed by Pearson and Bowmer"», the rate of 
growth of circulation being introduced by a phase lag, dependent on frequency 
parameter. However, the simpler approach adopted in this note has enabled the 
rate of growth of circulation to be included, while still avoiding computational 
difficulties. This does appear to result in a reasonable agreement between theory 
and experiment in the cases examined; also it throws up some additional qualitative 
points. For example, the physical difference between the “‘surge’’ of multi-stage 
compressors and the “‘surge”’ of single-stage compressors or pumps (i.e. compressors 
operating on an incompressible fluid as at low speeds) is not unknown experimen- 
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tally, although it has always been difficult to define formally. The surge of a multi- 
stage compressor at the speeds in question is often violent: it is certainly 
unmistakable. More particularly, the surge is usually manifest, if allowed to 
continue, by a series of thumps corresponding to a sudden pressure drop, a gradual 
recovery, followed by another drop, and so on. There are usually some } to 2 
cycles per second, often corresponding to the frequency of the whole system as a 
Helmholtz resonator. By contrast the surge of a single-stage compressor or pump 
is often difficult to determine, although this may be largely a matter of the pressure 
magnitudes involved. An essential difference is that the surge of a single-stage 
compressor or pump is characterised by a drop of pressure but no recovery. The 
surged flow is stable, although not usually axially symmetrical. 


One final comment: experience with part-load performance estimation by the 
method of Ref. 9 has shown that compressors are extremely fickle in the form of the 
curve of AT /(4N’) against N that they produce. Quite small changes in this curve 
could considerably affect its derivative and hence, if this hypothesis is correct, the 
shape of the surge line. Failure to obtain either qualitative or quantitative correla- 
tion of surge behaviour with the major performance parameters is not therefore 
surprising. 


The main object of this paper is a study of the effects of rate of growth of 
circulation on the stability criteria for axial flow compressors. Although supporting 
evidence has been produced, the theory cannot be regarded as verified. Hence no 
attempt is made to examine any practical consequences arising from the 
stability criterion. 


4. Conclusions 


A study has been made of the flow stability in multi-stage axial compressors. 
A hypothesis has been advanced from which a criterion for the stability of a stage 
has been established. The method of approach differs from previous work in that 
a finite rate of growth of circulation is postulated. It has also been found necessary 
to study the stability with respect to speed changes rather than the usual mass flow 
increment. The criterion so determined has been applied to multi-stage 
compressors. It has been concluded that instability will always originate in the last 
stage, even though that stage is unstalled. It has been necessary to introduce one 
empirical factor into the analysis but, when this is given a value compatible with 
other experimental and theoretical evidence, reasonably good agreement with experi- 
mental results is obtained. Additionally, some peculiarities noted in both the 
theoretical and experimental curves is strong prima facie evidence in favour of the 
hypothesis. Finally it is suggested that the surge of a multi-stage compressor is a 
different physical phenomenon from the surge or stall of low speed compressors 
or pumps. 
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Appendix I 


STABILITY CONDITIONS FOR A STAGE 


Four variables are used to define the operating point of a stage. These are the 
representative blade speed U, the fluid axial velocity V,, the pressure rise AP, and the 
throttle position which can be denoted by “t.” Of these quantities only two are 
independent. Let U and V, be the independent variables. Then for any small change 


where, in general, the throttle position would have to assume a value compatible 
with this change, given by a similar expression. 


Since 6U and 6V, are independent, they can be chosen in such a way that AP 
remains constant (i.e. 3AP=0). Equation (22) then becomes 


(57-), Gar), 


In an unstalled stage (@AP/O0V,)y is negative and (0V,/0U),p is positive; hence 
(dAP/dU), is always positive under these conditions. 


It can also be specified that 5U and 8V, be chosen in such a way that “t” is 
maintained constant. In this case equation (22) becomes 


= +GF), 
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Hence, substituting in (23), 


_ (OAP/AU)r, 


or. (@AP/aV.), (26) 
Hence if the stage remains unstalled the condition for stability 
The stage pressure rise AP will be given by 


where p is the density at inlet and » is the stage efficiency. In the unstalled range 1 
can be taken as constant, so that the condition (28) for stability becomes 


where the stage temperature rise is given by 


QU? (tan a, —tan a,) (31) 
gJK, (tan a, +tan 2,) 


“(tan a, —tan 


AE == 


glK, 
Differentiation with respect to U gives 


OAT, 2Q0U QU? (tan a,+tan a,) 0a, (32) 


Hence, considering this as taken at constant throttle and constant axial velocity con- 
ditions, the condition for stability can be written as 


provided that (tan a, +tan a,) > 0 Go 


The condition (34) will be satisfied in all normal axial compressors. 


Additionally, condition (33) may be deduced as an aerodynamic flow instability as 
follows. 
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With the stage operating in equilibrium suppose that the speed is suddenly increased 
by a small amount. Let the mean blade speed increase by 5U from U to (U+8U), 
Momentarily the mass flow, and therefore the axial velocity, will remain unchanged, 
The blade incidenc will consequently suddenly increase, and a new operating point 
not on the constant throttle line will be established. The change of incidence or fluid 
inlet angle will be given by 


a 
8U. 


As the mass flow, or axial velocity, increases to that demanded by the new speed, the 
inlet angle will decrease. It follows that 


5a, > Az, for stability . (36) 


where Az, is the ultimate change of fluid inlet angle to the constant throttle operating 
point. If this condition were not satisfied a state would be reached in which the actual 
operating fluid inlet angle was less than that for the equilibrium point. As the mass 
flow increased the inlet angle would decrease further, and equilibrium could never be 
achieved. 


We can write 


Then substituting in the inequality (36) gives 


ou /, 


0a, 


- (38) 


which is condition (33). 


If a finite rate of growth of circulation during a transient speed increase is assumed 
the condition (33) or (38) must be rewritten 


0a, 0a 
aU), for stability. (39) 


Now the circulation I’ round a blade is given by 


L’'=s(V,,, —Vw.)= sV, (tan a, —tan a,), (40) 
or 0a, 
whence sV, sec? a, +s (tan a, —tan z,) (41) 


From the velocity triangles at inlet 


U=V,tana,+V,tana, . . (42) 
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st) a, 4 
and so au Jy, (43) 
whence, substituting (4) with dV,=0, 
_ 
or or 
Now assume that Gu ps ‘ (45) 
Then equation (41) yields 
_ sec? a, (2,/0U), 46 
ou /, (tan a, —tan @,) (46) 
Also equation (42) gives 
ov, 1—V, sec? a, (da,/dU), 
oU ) (tan a, -F tan a,) (47) 
Hence, equating (46) and (47), 
2, —(1—p) tan a, +tan a, (48) 
tan a, +tan 2, mA 
and the condition (39) for stability becomes 
a, —(1 —p) tan 2, +tan a, cos? a, 


Appendix II 


STABILITY CONDITIONS FOR A MULTI-STAGE COMPRESSOR 
it will be assumed that the outlet throttle is choked, as at pressure ratios above about 
2:1. Then if all losses in the outlet ducting system are neglected the constant throttle 
line will be defined by 


Also W=pAV,, for continuity, . ‘ : (51) 
whence, for the last stage, Vi=k/ (Ty). (53) 
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At inlet to the last stage, 


U=V, (tana,+tana,), . ‘ (54) 


whence, considering the fiow in the ultimate state when all transients have disappeared, 
cos? a, 
= {1-7 su). } 


cos? a, 


(tan a, +tan 


where AT is the overall compressor temperature rise, gives for the last stage 


= {1- ar ( (57) 


Substitute in the condition for stability given at (45) of Appendix I; then 


{1- U {? tan z,—(1—p) tan a, +tan 
= a, +tan a, 
oAT a,+tan a, 


To evaluate the overall temperature rise at various speeds for constant throttle setting 
it is convenient to write 


=F+f(U) along the constant throttle line. . 
= 
A 
7 _ AT+3U*f (U) (62) 
7, Ty 
The condition for stability in (59) then becomes 
tan a,-+tan 1 AT+1U°fU 
0 1 meas (63) 
tana,+tana, ~(1—p) Ty 
4 Tn 
U stage AT +4U*f(U) (tan «,+tan 
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| Asymmetrically Plate-Reinforced Circular 


ed, | Hole in a Uniformly End-Loaded Flat Plate 
RAYMOND HICKS, Ph.D., M.Sc. 
(G.E.C. Simon-Carves Atomic Energy Group) 
(55) 


Summary: General expressions are obtained for the stress distribution in a flat 
56) plate end-loaded by equal principal stresses and containing a circular hole 
reinforced by an annular plate. The reinforcement is placed on one side of the 
plate only so that the other side remains flush. It is shown that the effect 
of an asymmetrically placed reinforcement is to induce, in the plate and rein- 
forcement, bending stresses which are not negligible compared with the 
externally applied direct stress. From particular examples it is found that 
57) reinforcements having practical dimensions can give rise to stress concentrations 
in the same order of magnitude as those in a similarly loaded plate containing 
an unreinforced hole. For practical purposes, data have been obtained which 
enable the designer to determine quickly the stress distributions associated 
with reinforcements and plates having a reasonably wide range of dimensions. 


58) 1. Introduction 


The problem of a symmetrically reinforced hole in an end-loaded plate has 
been considered by Timoshenko™, Gurney” and many others, for example 
in Refs. 3, 4, 5, 6. Timoshenko used a compact ring as a reinforcement, the 

_ radial width of which was small compared with the diameter of the hole. This 
ng | implies that, in the plane of the plate, the bending stiffness of the reinforcement is 
| small and the corresponding shear stresses negligible. Gurney used a relatively 

thin annular plate as a reinforcement and determined the stress distribution in 
| the reinforcement, as well as the plate, for various stress systems applied at infinity. 


0) 

Mansfield?) and the author‘ have considered the problem of a variably 

| reinforced circular hole in an end-loaded plate, the variation in the cross-sectional 

1) | area being chosen. so that the tensile resistance of the reinforcement increases at 
points around the hole where the stresses in the plate tend to be high. 


2) Elsewhere, the author’ has considered the problem of an asymmetrically 
ring-reinforced circular hole in an end-loaded infinite flat plate. That is, a 
reinforcement whose radial width is small compared with the diameter of the hole, 
so that it has no radial stress variation. In this case it is found that the asymmetry 

| of the reinforcement has two effects on the stress distribution in the plate. Firstly, 

3) bending stresses are induced in the plate and reinforcement. In general it is 

found that, for the plate, these bending stresses are less than the corresponding 

stresses in a plate with the relatively thin annular type of reinforcement considered 
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here. Secondly, for a given area of cross section, the asymmetrical reinforcement 
is less effective in reducing the stress concentration at the middle plane of the 
plate than a corresponding symmetrical reinforcement. This follows from the fact 
that elements of the asymmetrical reinforcement adjacent to the plate are elongated 
tangentially by bending of the reinforcement as well as by hoop stresses, and this 
results in less resistance to radial displacements of the plate around the hole. As 
shown later, the asymmetry of a plate reinforcement considered here has no effect 
on the stress distribution at the middle plane of the plate. 


This analysis deals with the problem of a circular hole in a flat plate uniformly 
loaded in its own plane by equal principal stresses applied at the boundary. The 
hole is reinforced by an annular plate having a radial width which is not small 
compared with the diameter of the hole. That is, there is a radial variation in the 
stresses acting at the middle plane of the reinforcement; furthermore, this variation 
is the same as the well-known radial variation for thick cylinders. 


The reinforcement is located on one side of the plate only (Fig. 1), so that 
the other side remains flush. This implies that bending moments are induced in 
the structure at the junction of the reinforcement and plate, since the middle plane 
of the plate does not coincide with the middle plane of the reinforcement. Also, 
since practical reinforcements may be several times as thick as the plate, it follows 
that the magnitude of these stresses may be high. That is, in the same order of 
magnitude as the stresses at the middle plane of the plate. 


As would be expected, for a plate under the combined action of bending and 
end load, the bending stress concentrations are not linearly related to the applied 
stress but become relatively more severe as the stress at the middle plane of the 
plate increases. For example, for a plate reinforced so that the stress distribution 
at the middle plane is the same as the stress distribution in a similarly loaded plate 
with no hole, the radial bending stress in the plate adjacent to the reinforcement is 
practically equal to the stress at the middle plane of the plate when the applied 
stress is 17,200 Ib./in*. That is, for this particular example, the resultant stress 
concentration in the plate with an asymmetrical hole is the same as the stress 
concentration in a similarly loaded plate with an unreinforced hole. 


A general solution is obtained for an annular plate with a fixed or simply- 
supported boundary. This solution is extended to deal with the problem of a 
similarly loaded infinite plate. It is found that, when the least lateral dimension 
of a plate is several timés the diameter of the hole, the solution for an infinite 
plate can be applied with sufficient accuracy to a plate of finite size. Furthermore, 
for any particular problem, the degree of approximation can be determined readily 
by comparing the stress applied to the boundary of the finite plate to the stresses 
acting around the corresponding contour in an infinite plate. 


For almost all practical problems numerical values for the stress concentrations 
can only be found by evaluating Bessel functions that are outside the range given 
in standard works of reference. It has been necessary, therefore, to tabulate data 
which will enable the designer to evaluate quickly the important stress concentrations 
for any particular problem within a reasonably wide practical range. 
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Figure 1. Section through plate and reinforcement. 


deflection of plate 

radial co-ordinate 

radial slope of plate 

thickness of plate 

inner and outer radii of annular reinforcement 
inner and outer radii of annular plate 
(thickness of reinforcement) / (thickness of plate) 
applied stress 

stress components at the middle plane of the plate 
moment components for the plate 

shear components for the plate 

Young’s modulus 

Poisson’s ratio 

flexural rigidity of plate 

constants defined by equation (2) 

constants defined by 2,h/D and 8,h/D 
dimensionless parameter r / « 

Bessel functions 

order of Bessel function, given by / (1+) 
arbitrary constants 

constants defined by equation (22) 

constants defined by equation (26) 

design parameters defined by equation (34) 
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FIGURE 2. Forces and moments interacting between plate and reinforcement. 


2. Theory 


The stresses acting at the middle plane of the plate (Figs. 1 and 2) are given 
by the following expressions 


B, 
r=a,+ 
pir (2) 
[p—(rr)o] 
For force equilibrium in a direction normal to the surface of the bent plate, 
,dw 


The negative sign has been inserted to satisfy the convention that downward 
deflections are positive. 


Using the classical theory for bending, equation (3) becomes 


d 1 dw 
d’w dw B,\, dw _ 


For convenience, let 


2=a,h/D,  B=B,h/D, 


x=rv 2, n’?=1+8. 
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Then equation (4) becomes 


This is a modified Bessel equation of order n; its complete solution is given by 
¥=C,I, (x) +C.,K,(x) (n integer) . ‘ (7a) 

For the range of problems considered here, n can be fractional or integer. However, 

it is proposed to use relation (7a) for all values of n, it being understood that when 


nis fractional K, (x) is given by 


sin nt 


Similarly, using the suffix ‘“‘R’’ to denote symbols associated with the reinforcement, 
itcan be shown that 


(rrp) = — (662) =2,2— Bin (9) 
r 
_  @b’ _ 2 
Also, for the reinforcement, 
Ve (Xr) (Xz), e (1 1) 


where n is fractional or an integer. 


3. Boundary Conditions 


3.1. GENERAL 
For compatibility of displacements between reinforcement and plate, 


Substituting for 66, etc., this gives 


(c? —b’) 
= 13) 
A\ Bb? -a — 
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Referring to Fig. 2, for moment equilibrium 


r 


that is, . 09 


Around the curve r=5, the radial slope of the plate is the same as the radial slope 
of the reinforcement, so that 


Also, since the radial moment is zero around the hole, 


The remaining boundary conditions depend on the manner of supporting the plate 
and are given in the following sections. 


3.2. SIMPLY-SUPPORTED PLATE 
For a simply-supported plate, 


dy 
=o. . « « 


In this case the boundary conditions (16), (17), (18), and (19) reduce to 


Dya{C,[ (6/2) + 2) Kv +5 (bv) |} - 
{Cr [ Tan! (BV an) + Inn (OV 20) | + 
+ Kan! (OV an) + Kum (BV 2x) | } + 
Cn (bV 2) + (bY 2) %2)=0 $ (20) 
Cay [ + a) | + 
+ [ Kas’ (a an) + Kan (@V/ ax) | =0 


(eva) + (eva) | +¢, [ Ke’ (eV 2)+ Ky (eva) | =0 


where the primes” denote differentiation with respect to x. 
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Solving these equations, it is found that 


1)TK 


Cu, = — A-1)SK 
TSK 
where 
ap) + Kran (a¥/ 29) | 
[1 (cy 2) + (ev) | 
[Ku (ev 9) + (ev ) | 
2) (by 2) 2) 


DV 2) |-S [Kv (bv 2) Ka (ov 2) 


{ [ (bV/ ap) + (bv 2x) | 


(6/20) + Kum 20) |} 


These expressions for the constants C, etc. are explicit, since the radial 
stress (rr), is defined by relation (13). It is interesting to note that, for the 
corresponding problem of an asymmetrical ring-type reinforcement (i.e. a reinforce- 
ment with a negligible bending stiffness in the plane of the plate), the radial stress 
(rr), cannot be obtained explicitly, since it is a function of the constants defining 
the deflection surface of the plate. Because of this the order of the Bessel function 
can only be found by solving a transcendental equation. This particular problem 
is discussed by the author in Ref. 9. 


Using the values just given for the constants, the radial and tangential moments 
in the plate are given by equations (23) which follow, where (rrp is given by 
telation (13). 
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[ky 
Mi=DVa(m K {[vh’@)+ - 


(23) 


-TS, [ (Xp) + (x) | } 


75x [ + (xe) ] 


3.3. Bur_t-IN SUPPORT 


For this particular problem the radial slope of the plate is zero at r=c, and 
it can be shown that expressions (21), (22), and (23) are valid provided that S is 
defined as 


(24) 


3.4. INFINITE PLATE 


The induced stress distribution in an infinite plate is defined by a modified 
Bessel function of the second kind. In general the absolute value of such a function 
decreases rapidly as its argument increases and is found to become negligible when 
the radial co-ordinate is several times its value at the junction of the reinforcement 
and plate. This implies that the solution for an infinite plate can be applied to a 
plate of finite size, provided that its least lateral dimension is several times the 
diameter of the hole. 


Proceeding as before, for an infinite plate, it can be shown that 


C,=0 
2DK, (bv 2) [LV aR] 
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Sz is defined in equations (22). 


r (26) 


(rr), = — 


4. Numerical Example 


To illustrate the order of the bending stresses induced in an infinite plate by 
the asymmetry of a reinforcement, consider the particular case where b=2a, 
a=30h and v=0-3. Also, for convenience of calculation, assume A=1-955. 
Then from (2), (5), and (13), it is found that 8=0, n=1, and (rr),=p. That is, 
the stress distribution at the middle plane of the plate is the same as the stress 
distribution in a similarly loaded plate with no hole. Also, for the particular 
parameters given here, numerical values for the Bessel functions are given in standard 
publications (for example, Refs. 10, 11, 12). 


For the reinforcement 


so that Ny” = 1 — (xg),”. (28) 
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Ficure 4. Example of tangential stress distribution. 


Assuming E =30 x 10° Ib./in.*, and for convenience taking p=17,200 Ib. /in.’, n is 
zero, so that the Bessel functions are given again numerically in standard 
publications. 


Thus, using (25) and (26), the stress distribution in the reinforcement and 
plate are as shown in Figs. 3 and 4. From these figures it is clear that the radial 
bending stress concentration in the plate, and the total tangential stress concen- 
trations in the reinforcement, are both substantially the same as the tangential 
stress distribution in a similarly loaded plate with no reinforcement, that is, 2p. 


5. Bessel Functions 


Except in a few special cases, it is found that, for practical reinforcements, the 
order of the Bessel functions and the magnitude of their arguments are outside 
the range of those given in standard works of reference. For this reason it is 
proposed to obtain numerical values for a set of parameters which can be used to 
determine the stress concentrations in plates and reinforcements having a reasonably 
wide range of practical dimensions. 

Considering the reinforcement, the order of the appropriate Bessel functions 
is given by 

n;? = 1 + 1- (6? —a¥ 
Thus, since the applied stress is positive, it follows that ng< 1. 


(29) 


Also, the order of the Bessel functions is real if 


(7) < « « 


12(1—v*) 
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Now the effect of increasing the reinforcement is to increase the ratio (rr)p |p, so 
that when the ratio is unity the stresses at the middle plane of the plate are the 
same as those in a similarly loaded plate with no hole. Since there is no practical 
advantage in increasing the reinforcement to give values of (rr)»/ p greater than 
unity, it may be assumed that, as an upper limit, (rr)y = pa. 


In this case, from (26) and (29), and assuming v=0°3, it is found that 


2 
(4) <7 07+ . « » 
R 


This relationship implies that, even for small values of E/p (say less than 2,000), 
the order of the Bessel functions is real, for practical values of a/b, provided that 


diameter of hole in reinforcement 
thickness of reinforcement 


< 30 approximately. 


In view of this result, design data will be obtained for all values of mg in the range 
0<ng<1. From (28), it follows that for the reinforcement the maximum value 
of Xp will be 


That is, xz is always sufficiently small to justify the use of the following series to 
determine the value of the Bessel functions. 


(4z)"t?" 


where v=+Np 
For the plate, x is not always sufficiently small to justify the use of the 


series (32). In this case values of K, (x) and K, (x) can be obtained from Ref. 10, 
while K, ,. (x) is given by 


1/2 


(=) 


Also, using the relationships 


(x) I, (x)] 


sin nt 


IG) =i-*J, (i). 


values of K,,, (x) and K,,; (x) can be calculated from information given in Ref. 11. 
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Finally, for intermediate values of n, values for the Bessel functions can be 
interpolated graphically. 


6. Design Data 


The stress distribution at the middle plane of the plate and reinforcement can 
be determined readily for any particular example, since (rr), is given explicitly 
by equation (13). Thus, for practical application, it is necessary only to tabulate 
data for the bending stress concentrations. In general the significant bending stress 
concentration in the plate acts in a radial direction and occurs at the point where 
the plate is attached to the reinforcement. Also, the significant bending stress 
concentration in the reinforcement occurs on the boundary of the hole and acts in 
a tangential direction. 


The design data in the rest of this section have been obtained so that these 


important bending stress concentrations can be evaluated readily for a wide range 
of practical plate and reinforcement dimensions. 


4:0- 
bel-3a 
b= 15a 
bs 2a 
O-4 av(«,) ie) 


FiGurE 5. ¢, for various values of av (a,) and b/a. 
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FIGURE 7. $4 for various values Of a (ag) and b/a. 
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For the plate and reinforcement, 


(Ms), = Dv 2C, | Ky’ (bV 2) + | 


= (av 25) |+ Cro [ (av a) + (33) 


Km(av en) |} 


For convenience in tabulating design data, the following parameters are introduced 
into the analysis: — 


[ Tax! 20) + [ (av an) + Kn 
(by 20) + Kan (by 20) + (av 20) | 


an) + =" — Tan (av 22) | (OV } 
K,, (bv 2) 


t= 


Using the appropriate series for the Bessel functions, numerical values for 
these parameters have been evaluated and are given graphically in Figs. 5, 6 and 7. 


With reference to Fig. 6, for (rr), / p less than one, Bessel functions defining 
¢p become imaginary for sufficiently large values of b/2. Thus, for b/2>3, 
points on the curves for ¢p have been obtained by using the approximate method 
now described. 


Figure 3 shows that the attenuation of the induced radial moment is very rapid 
in the plate, so that, when the radial co-ordinate is 1:55 (that is, 3a for this 
particular example), the magnitude of the bending stress is only 7 per cent of its 
value at the junction of the reinforcement and plate. 


In general, for any particular case, the attenuation of the radial stress becomes 
more rapid as the stress applied to the plate increases. Thus, for practical rein- 
forcements, and for sufficiently high values of the applied stress, the variation in 
the stress at the middle plane of the plate is small in the region of the plate where 
the radial bending stresses are concentrated. This implies that the local effects 
of the induced bending can be determined by assuming that the plate is everywhere 
under the action of constant radial and tangential stresses of magnitude (rr)». 
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In this case B=0, n=1, and 


K, () 
[x ‘ ‘ ‘ (35) 
7 | 


where y=b [2 


For any particular example, the accuracy of this simplification increases as by 2 
increases and any error will result in an over-estimate of ¢p. 


From equations (25), (26) and (34) 


K, (6/2) L 

or= 
by an | Ky (bv 2) + Ka (bv) |P 
Lia 
dp Pv ap 
ap 
PV Ly a—A ag 

60=P | 25) + Kus (av 


Thus, in terms of these parameters, the moment concentrations become 


_ Groh? (A —1) on 
(M,), 2 (or Ap) 


(A 1) (1 v2) op 
2 (or —A° or) $0 


(37) 


For any particular example these can be evaluated readily from the curves given 
in Figs. 5, 6 and 7. 
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